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ABSTRACT 



Aims. To establish the chemical origin of a set of complex organic molecules thought to be produced by grain surface chemistry in 
high mass young stellar objects (YSOs). 

Methods. A partial submillimeter line-survey was performed toward 7 high-mass YSOs aimed at detecting H2CO, CH3OH, CH2CO, 
CH3CHO, C2H5OH, HCOOH, HNCO and NHjCHO. In addition, lines of CH3CN, C2H5CN, CH3CCH, HCOOCH3, and CH3OCH3 
were observed. Rotation temperatures and beam-averaged column densities are determined. To correct for beam dilution and deter- 
mine abundances for hot gas, the radius and H2 column densities of gas at temperatures >100 K are computed using 850 //m dust 
continuum data and source luminosity. 

Results. Based on their rotation diagrams, molecules can be classified as either cold (<100 K) or hot (>100 K). This implies that 
complex organics are present in at least two distinct regions. Furthermore, the abundances of the hot oxygen-bearing species are cor- 
related, as are those of HNCO and NH2CHO. This is suggestive of chemical relationships within, but not between, those two groups 
of molecules. 

Conclusions. The most likely explanation for the observed correlations of the various hot molecules is that they are "first generation" 
species that originate from solid-state chemistry. This includes HjCO, CH3OH, C2H5OH, HCOOCH3, CH3OCH3, HNCO, NH2CHO, 
and possibly CH3CN, and C2H5CN. The correlations between sources implies very similar conditions during their formation or very 
similar doses of energetic processing. Cold species such as CH2CO, CH3CHO, and HCOOH, some of which are seen as ices along 
the same lines of sight, are probably formed in the solid state as well, but appear to be destroyed at higher temperatures. A low 
level of non-thermal desorption by cosmic rays can explain their low rotation temperatures and relatively low abundances in the gas 
phase compared to the solid state. The CH3CCH abundances can be fully explained by low temperature gas phase chemistry. No cold 
N-containing molecules are found. 

Key words. Astrochemistry - Line: identification - Methods: observational - Stars: formation - ISM: abundances - ISM: molecules 



1. Introduction 

A characteristic stage of high-mass star formation is the hot core 
phase where high abundances of complex molecules are found in 
the gas-phase and in the ice mant les around the dust grain s in the 
cold extended envelope (e.g. iMillaii 119961: llkeda et all I2OOII: 
' iBoogert & Ehrenfreundl l2004t) . The existence of both ices and 
gas-phase molecules around the same high-mass protostars sug- 
gests that the envelopes have temperatures that range from >200 
K in the inner region to as low as 10 K in the co l d outer part 
(e^ lvan derTak etali 120001: iGibb et all l2000bt llkeda et all 
l200lh . Together the submillimeter and infrared data indicate 
that there is a rich chemistry which involves interactions be- 
tween gas and grain-surface species. The precise origin of com- 
plex molecules is still heavily debated: are they "first gener- 
ation" spec ies formed on t he grains during the cold collapse 
phase (e.g.. lTielens & Charn lev, 1997), or "second generation" 
species produced by high-te mperature gas phase reactions fol- 
lowin g ice evap oration (e.g.. lMillar et aUri991l : ICharnlev et al.L 
1 1 9921: ICharnlevl 1 1 995h ? The aim of this work is to test the for- 
mation schemes of key complex organic molecules by presenting 
submillimeter observations of a sample of high-mass cores using 
the same instrument and frequency settings. 

High column densities of the organic molecules CH3OH, 
CH4, HCOOH, OCN" and potentially H2CO have been de- 
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tected in ice mantles in the cold outer envelopes of high-mass 
young stellar objects (YSOs) (Gibb et al., 2004). Abundances 
for some ice species such as CH3OH and OCN" are large, 
up to 30% with respect to H2O ice or 3x10"^ with respect to 
H2 (B rooke et al., 199 9:lDartois et al LI 19991: iPontoppidan et all 
\200?i Ivan Broekhuizen et al.L l20'ol . Information on minor 
species such as CH4, HCOOH, H2CO and CH3CHO is more lim- 
ited, but they appear to have co nstant abundances at the level of a 
few % with respect to H 2O ice dSchutte et aUll999l:lKeane et all 
l200ll : lGibb et al.Ll2004l) . More complex species with abundances 
lower than 0. 1% cannot be detected by infrared observations be- 
cause of their small optical depth. Moreover, features of the same 
functional groups such as C-H and O-H groups overlap, so that 
they can only be assigned to a class of molecules rather than to 
a specific species. Deep submillimeter surveys in the gas phase 
are more sensitive to detect these minor species. 

The abundances of the detected ice species are so large that 
they cannot be produced by pure gas-phase processes but must 
be formed on the grains. A proposed scheme for the forma- 
tion of various "first gener ation" molecules on grains (based on 
[Tielens & Charnlevl 119971) is shown in Fig.[T] Through succes- 
sive addition of H, C, O, and N to CO, many complex organic 
species are formed. The outcome depends on the relative gas- 
phase abundances of the atoms, their mobility on the grains, and 
the CO abundance. If the formation routes depicted in Fig. [1] are 
indeed valid and proceed under similar initial conditions, this 
would imply that the molecules in this scheme are related and 
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that they will have correlated abundance ratios upon evapora- 
tion. For different initial conditions, the relative abundances will 
vary from source to source; the H2CO/CH3OH ratio depends for 
example on ice temperature, and the CaHgOH/CHsOH ratio on 
the atomic carbon abundance. In contrast, gas phase formation of 
second generation species is expected to cause time-dependent 
abundance variations where "first" and "second generation" as 
well as non-related "first generation" species are uncorrected. 
A combination of rotation temperatures and relative abundances 
between species can determine which species are part of the 
same chemical group. To test the basic grain surface chemistry 
scheme, our survey has targeted emission lines of all species de- 
picted in Fig.[T] Only CO2, the unstable species HCOO, HCCO, 
and NH2CO, as well as CH2(OH)2 and NH2CH2OH for which 
no laboratory frequencies are available, are lacking. 

Many line surveys have been performed to study the molec- 
ular composition of high mass star forming region s (e.g., 
Blake et al.Lll987l:ISutton et alill995l:lHelmich & van Dishoeck , 



1997HSchilke et al.L 119971: iNummelin et all l2000t iGibb et al 
2000a ISchilke et al.Ll200l . However, they often cover a sin 
gle atmospheric window, each at a different telescope. These 
data thus probe a variety of excitation conditions with differ- 
ent beam-sizes, which complicates the comparison of column 
densities and abundances between different sources. Our observ- 
ing strategy employs deep, partial line surveys in multiple at- 
mospheric windows rather than complete, unbiased line surveys. 
Partial line surveys exist, but are usually targeted at a few spe- 
cific m olecules such as CH3O H jvan der Tak et al.L I2000D and 
HNCO (IZinchenko et al.L [2000h . The advantage of our survey is 
that it efficiently covers lines with a large range of excitation 
temperatures and that the settings and telescopes are the same 
for all sources. The lines of the molecules targeted in this work 
are sufficiently strong so that no complete line survey is needed 
for their iden tification. Compared with previous studies such 
as those by Ikeda et alj (HOOI) and Zinchenko et alj (12000') our 
observations cover higher frequencies, which implies smaller 
beams and, for some species, lines from higher energy levels. 

This paper is organized as follows: Sect. |2] presents the se- 
lection of sources plus frequencies, and the observational tech- 
niques, Sect.[3]presents the data reduction techniques, Sect.|4]the 
results on the derived rotation temperatures plus column densi- 
ties and compares those to other high mass YSOs. In Sect. |5] 
the correlations of the abundances between different molecu- 
lar species are determined. Sect. |6] discusses the origin of the 
constant abundance ratios and compares these to astrochemical 
models, and finally in Sect. [7] the main conclusions are summa- 
rized. 
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Fig. 1. Grain-surface chemistry routes involving hydrogenation 
of CO. Solid rectangular boxes contain molecules which have 
been detected in interstellar ices, whereas dashed boxes indicate 
molecules that have been detected in the gas phase (based on 
lTielens&Charnlevlll997h . 

Table 1. Positions, luminosities, distances and typical line width 
of selected sources. 



Sources 


Q'(2000) 


(5(2000) 


L" 


d' 


AV 








[10= Lo] 


[kpc] 


km s"' 


AFGL2591 


20:29:24.6 


+40:11:19 


0.2 


1.0 


3 


G24.78 


18:36:12.6 


-07:12:11 


7.9 


7.7 


6 


G75.78 


20:21:44.1 


+37:26:40 


1.9 


4.1 


4 


NGC 6334 IRS 1 


17:20:53.0 


-35:47:02 


1.1 


1.7 


5 


NGC 7538 IRS 1 


23:13:45.4 


+61:28:12 


1.3 


2.8 


4 


W 3(H,0) 


02:27:04.6 


+61:52:26 


0.2 


2.0 


6 


W33A 


18:14:38.9 


-17:52:04 


1.0 


4.0 


5 



, except for G24.78 and G75.78 where 
1 ( l2005h and lChurchwell et alJ d 19901) 



" From Ivan der Tak et al] ( 12000 
the luminosities are from|Zhang 

respect ively. 

From lvan der Tak et all (|2003), exce pt fo r G24.78 and G75.78 wh ich 
are based on 
spectively. 



Forster&Caswelll ( ll989h and lChurchwell et alj jl990h re- 



2. Observations 

2. 1 . Source sample 

The basis for this survey is formed by a sample of 7 high-mass 
young stellar objects summarized in Table [T] The sources are 
known to have hot gas in their inner envelope where ices have 
evaporated and they either show complex organic molecules 
due to hot core chemistry or a re on the verge of reaching that 
stage (e.g.. lBoonman et 31112001 ). Most of the sources were se- 
lec ted from the sample of high mass YSOs studied in detail 
by Ivan der Tak et al.l (l2000h based on the following criteria: (i) 
narrow line widths (<6 km s"') to minimize line confusion; 
(ii) detection of CH3OH lines, where the strongest lines in the 
7a:-6a' 338 GHz band should have main-beam temperatures of 
>1 K; (iii) distances <5 kpc and (iv) visible from the James 



Clerk Maxwell Telescope (JCMTfl Preference was also given to 
sources for which infrared absorption line data of ices are avail- 
able. The source list additionally includes G75 .78 and G24.78 . 
G75.78 has been previously observed by Hatch ell et alj ( |1998[) 
and obeys all previously mentioned selection criteria. G24.78 
has a distance >5 kpc, but was sele cted based on the large line 
intensities of CH3CN detected by iBeltran et all (HOOS). Trial 
JCMT observations of the 7k-6k 338 GHz CH3OH band con- 
firmed that the strongest lines are > 1 K and G24.78 was there- 
fore added to the sample. 



' The James Clerk Maxwell Telescope is operated by the Joint 
Astronomy Centre, on behalf of the Particle Physics and Astronomy 
Research Council of the United Kingdom, the Netherlands Organization 
for Scientific Research and the National Research Council of Canada. 
Project IDs are m04bnl2, m05an05, and m05bnll. 
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2.2. Observed frequency settings 

The observed frequencies are shown in Table |2] and were se- 
lected based on whether they contained at least one strong line 
for a target molecule as well as lines of other interesting species. 
Strong lines of target molecules were chosen because of their 
high main-beam temperat ures and minimum line confusion in 
line surveys of Orio n A bv lSutton et alj (1 19851) and Orion KL by 
ISchilke et all (llQg? *) at 230 GHz and 345 GHz respectively. The 
aim was to cover at least two transitions for each molecule to 
determine rotation temperatures: one with >100 K, and one 
with <100 K. Table|2]shows the additional lines present in the 
rest frequency side band. The molecular species detected in the 
image band differed, because it varied per source. Finally, for 
H2CO, CH3OH, and HNCO lines of the isotopic species were 
observed to determine the optical depth of the main isotopo- 
logues. 

An additional number of settings have been observed 
for CH2CO, CH3CHO, and HCOOH with the Institut de 
Radioastronomie Millimetriqu43 (IRAM) 30 m telescope for 
AFGL 2591, G75.78, NGC 7538 IRSl, and W 3(H20). The 
settings were chosen because of their high line strengths and 
low excitation temperatures. Combined with the JCMT settings, 
they allow for the determination of any cold component of these 
molecules. 

2.3. JCMT observations 

The frequencies given in the first part of Table |2] were ob- 
served for the 7 sources with the JCMT on Mauna Kea, Hawaii 
from August 2004 to January 2006. The beam size (Oheam) is 
20-21" for observations in the 230 GHz band and 14" in the 
345 GHz band. The spectra were scaled from the observed an- 
tenna temperature scale, T^, to main beam temperatures, Tmb, 
using main beam efficiencies /7MB 

of 0.69 and 0.63 at 230 GHz 
and 345 GHz, respectively. The front-ends consisted of the 
facility receivers A3 and B3; the back-end was the Digital 
Autocorrelation Spectrometer (DAS), covering 500 MHz in- 
stantaneous bandwidth with a spectral resolution of 312.5 and 
625 kHz respectively. Pointing was checked every 2 hrs or when- 
ever a new source was observed and was always within 3". 
To subtract atmospheric and instrumental backgrounds beam 
switching with a chop throw of 180" was used. We aimed for 
a Tims ~20 mK on the scale when binned to a ~1 km s ' 
channel, which was obtained with integration times of ~1 hr 
at 230 GHz and 2 hrs at 345 GHz dependent on weather For 
NGC 6334 IRSl the confusion limit was reached after 30 min- 
utes of integration, however, at which point the observation was 
terminated. The absolute calibration is accurate to better than 
15% from comparison to spectral line standards. 

The JCMT B3 345 GHz receiver was operated in single 
side band mode to minimize line confusion. The A3 230 GHz 
receiver does not have the option of i mage sideband supp res- 
sion. For line-rich sources like Orion dSutton et al.L Il985h . it 
is difficult to assign a line to the correct sideband. For our 
sources the 230 GHz spectra are not confusion limited except 
for NGC 6334 IRSl (see Fig. |3]l. For NGC 6334 IRSl, pairs of 
spectra with 10 MHz offsets were taken. This shifts lines in the 
image sideband by 20 MHz with respect to those observed in the 
rest sideband making a unique identification possible. 



^ The IRAM 30 m telescope is operated by the Institut de 
Radioastronomie Millimetrique on behalf of the Centre National de 
la Recherche Scientifique of France, the Max Planck Gesellschaft of 
Germany, and Instituto Geografico Nacional of Spain. 



2.4. IRAM Observations 

The frequencies given in the second part of Table |2] have been 
observed with the IRAM 30 m telescope. The beam size (flbeam 
in arcsec) is 29" at 86 GHz, 17" at 140 GHz, and 10.5" at 
235 GHz. The front-ends were the BlOO, CI 50, B230, and C270 
facility receivers. The VESR4 auto-coiTelator was used as the 
back-end, with a bandwidth of 120 MHz, and a channel spacing 
of 80 kHz. Pointing was checked every 1.5 hrs and was always 
within 2-3". Atmospheric and instrumental backgrounds were 
subtracted through "wobbler" switching where the telescope al- 
ternately observes positions +220" and -220" away from the 
source. Integration times differed per source and depended on 
the atmospheric conditions, but were generally 30 minutes for 
AFGL 2591 and G75.78, and 1 hr for NGC 7538 IRSl and 
W 3(H20). All observations were converted to the main-beam 
temperature scale using Fetf/Seff of 0.78 at 86 GHz to 0.46 
at 260 GHz. Conveniently, all observations were performed in 
single-sideband mode. 



3. Data analysis 

3.1. Rotation diagrams 

Line assignments were made by comparing the observed fre- 
quencies corrected for source velocity with the JPL0, CDMS0 
and NIS10 catalogs. Possible assignments had to be within 
~1 MHz, unless the uncertainty of the frequencies in the catalogs 
was larger For NH2CHO, for example, lines were only assigned 
if multiple transitions were detected in the same spectrum. 

Integrated intensities for the identified lines are given in 
Tables IB.H4B.15l These were used to derive rotation tempera- 
tures and column densities th rough the rotation diagram method 
jGoldsmith & Langeii [19991) when three or more lines were de- 
tected over a sufficiently large energy range. The integrated 
main-beam temperatures are then related to the column density 
in the upper energy level by: 



_ 3^/ TuBdV 



(1) 



where A^u is the column density of the upper level , is the 
degeneracy in the upper level, k is Boltzmann's constant, v is the 
transition frequency, fi is the dipole moment, and S is the line 
strength. The total beam-averaged column density A^t in cm"^ 
can then be computed from: 



Nt 



-EJT,„ 



gu Q(Tmt) 



(2) 



where 2(7^iot) is the rotational partition function, and is the 
upper level energy in K. 

The calculated rotation temperatures and column densities 
can be used to make predictions for the intensities of other lines 
for a specific species. This was used as an independent check 
on line-assignments. If lines predicted to have high intensities in 
the observed frequency ranges are not seen, a detection is con- 
sidered questionable. In the analysis, a coiTection has been made 
for the different beam-sizes at different frequencies. The assump- 
tion has been made that the maximum source size is equal to the 
smallest observed beam-size of the JCMT for all sources, i.e. 
14". The calculations have also been performed assuming that 



^ http : // spec . jpl . nasa . g ov/ftp/pub/cat alog/catform . html | 

|http : //www . ph i .uni-ko eln . de/vorhersagen/] 
' |http : //physics . nist . gov/PhysRefData/Micro/Html/ contents . html | 
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Table 2. List of observed frequency settings and molecular lines covered at the JCMT and IRAM. 



Molecule 


Freq. 






Transition 


Additional molecules 




[GHz] 


[K] 


[D2] 






JCMT 


HCO 


346.725 


41.61 


7.341 


4q 4 5 4— 3o,3,4,3 


C2H5CN/C2H50H/H2C"'0/HCOOCH3 




348.778 


172.36 


5.313 


42,3,4,4— 32,2,3,3 


HCOOCH3/CH3'^CN/C2H50H 


H2CO 


218.222 


20.96 


16.2985 


3o,3— 2o,2 


CH30H/C2H5CN/nCOOCH3 




364.275 


158.42 


52.1647 


53,3- 43,2 


C2H,OH 


CH3OH 


338.344 


5.3431 


69.36 


73,4,0— 63,4,0 


HCOOCH3 


CH2CO 


240.187 


88.01 


71.489 


12], 12— 1 1 1,1 1 


CH3 OH/H2C ' O/C2H5 CN/CH3 CN/HCOOCH3 




220.178 


76.48 


65.434 


lliai- lOuo 


CH3 ' 3 CN/" CH3 OH/HCOOCH3 




222.298 


182.26 


61.0830 


113,9- 103,8 


CH3 OCH3/CH3 CCH/TICOOCHi 


CH3CHO 


226.590 


71.30 


77.779 


12o,i2— llo,ii 


HCOOCH3 


C2H5OH 


230.991 


85.53 


48.953 


165,11- 164,12 


CH30H/C2H)5CN 




347.916 


250.88 


66.314 


2O4.17.0- 194.16.0 


C2H,CN/HCOOCH3/H2C'*'0 




339.979 


57.89 


17.724 


94,6- 83,5 


C2H5CN 


HNCO 


219.798 


58.02 


28.112 


10o,10,ll- 9o,9,10 


H2''CO/C2H5CN 




352.898 


187.25 


43.387 


I61, 15,17- 15l, 14,16 


C2H5CN/C2H5OH/HCOOCH3 


HN"CO 


240.881 


112.53 


30.4313 


111, 11, 12- lOl, 10,11 


CH3 OCH3/CH3 GH/TINCQ 


NH2CHO 


339.904 


245.79 


539.342 


166,11- 156,10 


same setting as C2H5OH 




340.491 


165.59 


605.520 


163,14- 153,13 


same setting as CH2CO 




345.183 


151.59 


664.219 


17o,17- I60.I6 


HCOOCH3/C2H50H/'3CH30H 


HCOOH 


225.238 


88.02 


17.734 


103,8-93,7 


H2CO/HCOOCH3/CH3 OCH3/" CH3 OH 




338.144 


180.53 


27.1548 


154,12- 144.11 


CH3OH/C2H5OH 




346.719 


143.92 


28.726 


152,13- 142,12 


same setting as HCO 




356.137 


158.67 


30.6650 


I62.15- 152,14 


NH2CHO/HCOOCH3/C2H5 CN/H2 ' ^ CO/CH3 OH 


IRAM 


CH2CO 


81.586 


22.86 


22.4974 


4l,3-3l,2 






140.127 


39.97 


41.1295 


7l,7- 61,6 






244.712 


89.42 


71.4841 


12i,ii- 111. 10 






262.619 


140.63 


25.3788 


132,12- 122.11 




CH3CHO 


98.901 


16.51 


31.2087 


5l,4-4l,3 






112.249 


21.14 


37.9261 


61,6- 5 1,5 






149.507 


34.59 


51.1868 


81,8- 7l,7 






168.093 


42.66 


57.7848 


9l,9- 81,8 




HCOOH 


223.916 


71.93 


18.7011 


102,9-92,8 






247.514 


150.71 


17.0080 


1 15,6/7- 105,5/6 






257.975 


83.91 


23.2080 


12l,12- 111, 11 






262.103 


82.77 


23.2998 


12o,i2- llo,ii 





the emission is extended with respect to the largest JCMT beam- 
size of 21". In this case, the rotation temperatures are higher 
and have larger errors by a factor of 2. In reality, as will be dis- 
cussed in Sect. 13.21 the hot core gas will be even more severely 
beam-diluted. Finally, the calibration uncertainties are taken into 
account in the rotation diagram analysis and are shown as error- 
bars on the rotation diagram plots. 

When only very few lines were detected or only a small 
range in was covered, no rotation temperatures or column 
densities could be derived. If the detection of the species was 
considered real, the column density was estimated assuming the 
rotation temperature to be the average for the same molecule in 
other sources. If no lines were detected, a 3cr upper limit was 
determined. For lines from isotopologues, rotation temperatures 
equal to the main species were used to estimate the column den- 
sities. 

A number of assumptions are implicitly made in the rota- 
tion diagram method. In short, they are that the excitation can 
be characterized by a single temperature Tiot and that the lines 
are optically thin. Even though the first assumption will almost 
certainly be incorrect, the rotation diagram method does give the 
average excitation temperature of the region from which most 
of the molecular emission arises. If the lines are sub-thermally 
excited because the density is below the critical density (nci), the 
rotation temperature is a lower limit. This could also be partly re- 



sponsible for the scatter in the rotation diagrams. Since for most 
molecules studied here no collisional rate coefficients are known, 
a more sophisticated statistical equilibrium analysis is not pos- 
sible. The critical density, n^, for the CH3OH 1(,.k-(!>6,k transi- 
tions is ~10^ cm"-' while for the H2CO 3o,3-2().2 transition it is 
~5xlO^ cm"-'. Since «„ is proportional to fi^ and v^, transitions 
with high values for either will be most affected. Most species 
in this study have dipole moments fi close to the CH3OH value 
of 1.7 Debye. H2CO, CH3CHO, HNCO, NH2CHO, and CH3CN 
have dipole moments 2-3 times larger and for these species sub- 
thermal effects should thus be more important. Optical depth 
can significantly influence the results as well. Column densities 
derived from optically thick lines will be severely under- and 
rotation temperatures overestimated. Lines of isotopic species 
can be used to determine the optical depth. Alternatively, it is 
possible to predict which lines are optically thick based on the 
rotation diagram results and exclude those from the fit. In this 
paper both approaches have been used. A final uncertainty is 
infrared pumping. It is a well-studied phenomenon for HNCO 
(IChurchwell et aI.LIl986 ). Its importance for other species is un- 
known, although all molecules have dipole-allowed mid-infrared 
transitions. However, no vibrationally excited lines have been 
detected in this survey except for CH3OH. With these limitations 
the rotation diagram method gives a useful indication of the tem- 
perature region where the abundance of a molecule peaks. 
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3.2. Beam dilution correction and iiydrogen column density 
in hot gas 

If the emission from the complex organic species comes from a 
warm region where the temperature is higher than the ice evap- 
oration temperature of 100 K, it becomes important to properly 
take the amount of beam dilution into account, especially for 
inter-comparison of the sources. The region where the tempera- 
ture is above 100 K is expected to be small (~ 1") compared to 
the size of the single-dish beam, and is related to the luminos- 
ity L. Also, the beam-dilution depends on the distance as d^. To 
correct for the amount of beam dilution we established a set of 
models for the dust envelopes of our sources based on SCUBA 
850 fim datcQ. For simplification we fix the inner radii of the 
envelopes to 500 AU, the outer radii to 100,000 AU (~ 0.5 pc) 
and adopt a power-law density profile, n oc r"'^, typical of a 
free-falling core. is chosen to be so large that gas extending 
beyond this radius will not significantly contribute. We assume 
that the total luminosity quoted in Table [1] is provided by a sin- 
gle 30,000 K blackbody at the center of the envelope and fit the 
observed 850 /urn fluxes of each source by adjusting the mass 
of the envelo pe using the DUSTY dust radiative transfer code 
dlvezic et all 1 1999) as described in |j0rgensen et alj ( l2002h . The 
results of these models provide a temperature profile for each 
source from which the radius where the temperature increases 
above 100 K (^7=100 k)^ and the beam-dilution can be estimated 
(see TableO. Figure|2]compares the inferred source sizes of the 
warm gas in arcsec with luminosity of the sources with all dis- 
tances set to 1 kpc. These correlate very strongly with luminos- 
ity as expected. Our results are consistent with the independently 
derived relation from dust modeling of a large range of sources, 
■'^r=ioo K ~ 2.3 X 10'"*( V-L/i©) cm (Doty private communica- 
tion). In fact, the luminosity itself can be used to estimate the 
linear size of the region where the temperature is higher than 
100 K. We therefore adopt these radii for the beam dilution cor- 
rection for molecules that unambiguously probe warm > 100 K 
gas estimated from the rotation diagrams. The derived beam di- 
lution (?7bf). 



(^r=ioo K/d)^ 



(0beam/2)2 + (Rt=IQO d d)^ 



(3) 



is then used to calculate source-averaged column densities 
(A^Taioo k) of hot molecules (X) in Sect. 14.41 



T>1{)() 



k(^) 



A^t(^) 

?7bf 



(4) 



The radiative transfer models used to calibrate the ^7-=iook - 
L relation can also be used to calculate the hydrogen mass in 
the inner warm regions and thereby provide an H2 column den- 
sity to determine molecular abundances for T > 100 K. With 
the derived density and temperature profiles listed in Table |3] 
the number of H2 molecules, MoZ(H2), within ^r=ioo k, can be 
calculated from: 



Mo/(H2) 



A-nr^rii \ — 



in gas with T > 100 K, Nn 
Mo/(H2) divided by nR^^^^^ , 



100 k( H2) in cm ^ is then simply 



Nt>ioo k(H2) - 



4nir" 



(3 - a)Rl^ 



(^7-=100 



K 'i J 



(6) 



100 K 



Since both Nt>ioo k(H2) and Nt>ioo k(X) refer to column densi- 
ties for the same temperature region, they can be used directly to 
calculate abundances for a species X in the hot gas through: 



[X] 



r>ioo 



k(X) 



Nt>\00 k(H2) 



(7) 



It is clear from Table [3] that the derived values of Rt=ioo k 
and Nt>ioo k(H2) are rather similar for most sources with the 
exception of AFGL 2591 and G24.78. AFGL 2591 has a signif- 
icantly smaller and G24.78 a much larger radius and H2 column 
density. These are related to their low and very high luminosi- 
ties. The beam-dilution for G24.78 is thus very similar to that 
for the other sources in our sample. This makes a good compar- 
ison between G24.78 and the other sources possible. 

Some of our sources have previously been modeled in 
more detail (e.g., |y an der Ta k et al 1 120001: iMueller et all I200I 
iHatchell & van der Tab, i2003l) : these models show some diff'er- 
ences, e.g., in the slope of the density distribution from source 
to source and e ven for the same s ource depending on method. 
As an example, Ivan der Tak et al. I (I2000h found that Rt=iqo k is 
5.3x10^ AU forNGC 6334 IRSl and 2.8x10^ AU for W 3(H20) 
with a density slope a-2. Both values are within 10-15% of our 
model result. Together with the good agreement with the models 
by Doty (private communication) this implies that the model- 
dependent uncertainty has a rather small effect on Rt=\oo k and 
the column densities. We therefore consider our adopted ap- 
proach to provide the internally most consistent set of estimates 
- even though we caution that a full self-consistent set of models 
(explaining, e.g., both single-dish and interferometric continuum 
data) is required to estimate abundances to better than a factor 2- 
3. 

Finally, interferometric observations have shown that even 
on scales as small as a few thousand AU, there may be 
strong chemical diff'erentiation, e.g., oxygen-bearing molecules 
peaking at a diff erent position than n itrogen-bearing s pecies 
(Liu, 2005; Blakeilill fT996l fT987t IWvrowski et all 119991: 
fSutton et al.L 119951) . From single-dish observations of a single 
source, it is not possible to determine whether the emission of 
various species arises from different regions. A comparison of 
the abundances of the species studied in this paper for a num- 
ber of sources does, however, show whether the emission likely 
comes from the same or a different region (Sect.|5]l. 



where ni is the density of H2 molecules at the inner radius, and 
a the power of the density profiles. The column density of H2 



* The SCUBA 850 /jm data were retrieved through the JCMT archive 
as a Guest User at the Canadian Astronomy Data Center, which is op- 
erated by the Dominion Astrophysical Observatory for the National 
Research Council of Canada's Herzberg Institute of Astrophysics 



4. Results 

4.1. General comparison between sources 

Figure [3] shows the spectra at 240.25 and 353.05 GHz toward all 
sources and Tables IB.lHB.15l contain all measured line intensi- 
(5|ies for the molecules discussed in this paper. There is a clear 
difference in the "richness" of the sources in molecular lines. 
The intensities and thus also column densities of "line-poorer" 
sources are lower for most organic species for sources such as 
AFGL 2591, G75.78, NGC 7538 IRSl, and W 33 A. Some lines, 
such as the HCOOCH3 314,27- 304,26 at 352.922 GHz, are weak 
relative to the HNCO 161.15- 15i,i4 at 352.898 GHz for the 
"line-poorer" sources, but strong for the "line-richer" sources, 
i.e. G24.78, NGC 6334 IRSl, NGC 7538 IRSl, and W 3(H20). 
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Fig. 2. Comparison of luminosity with source size where 
T =100 K. Rt=\()q k is given in arcsec assuming all sources to 
have a distance of 1 kpc. The solid line represent the best fit to 
the data. 

Table 3. Parameters of temperature and density profiles, and 
beam-filling factors derived from 850 yum SCUBA data and 
source luminosities. 



Sources 



AFGL 2591 

G24.78 

G75.78 

NGC 6334 IRSl 
NGC 7538 IRSl 
W3(H,0) 
W33A 



n" ^r=iooK 

,-3l 



[cm 



[AU] 



Nthoq k(H2) 
[cm-2] 



8.5(6) 
1.0(8) 
2.1(7) 
3.8(7) 
3.4(7) 
6.2(7) 
4.1(7) 



1.8(3) 
1.3(4) 
5.6(3) 
4.7(3) 
4.9(3) 
2.4(3) 
4.5(3) 



0.40 

4.8 

0.98 

1.8 

1.6 

2.9 

1.9 



7.6(22) 
4.0(23) 
1.2(23) 
2.4(23) 
2.1(23) 
1.8(23) 
2.6(23) 



6.2C 
5.8C 
3.70 
1.30 
6.0( 
2.90 
2.60 



" Density at inner radius of 500 AU. 
* The optical depth t is given for the dust at 100 fim. 
T/bf is the beam-filling factor. 



There therefore seems to be an intrinsic chemical difference be- 
tween the "line-rich" and "line-poor" sources. 

The two settings shown in Fig. [3] were aimed at detect- 
ing the CH2CO 12i,i2-lli,ii transition at 240.187 GHz and the 
C2H5OH 211,20,2-262,19,2 line at 352.858 GHz. Both transitions 
were detected in most sources, indicating that the selected set- 
tings are suitable for obtaining detections or otherwise limits. 
Many lines of other target molecules have been detected in these 
settings as well. 



4.2. Optical depth determinations 

Some species such as CH3OH, H2CO, and CH3CN have previ- 
ously been_joimdjo_b£_o£tical in their main lines (see 
e.g., iHelmich & van Dishoecki Il997h . Lines for selected iso- 
topic species were used to determine whether the lines of the 
main species are optically thick and to derive column densities 
for the main isotopologues (see Tables IB. 11 IB. 21 IB. 31 IB. 81 and 
IB. 10b . For H2CO and HNCO, the same transition could be used 
for different isotopologues. For CH3CN, the optical depth was 
determined by comparing the column density derived from the 
rotation diagram analysis for the main species and that derived 
from the detected transition for the isotopologues. The ratios for 
the beam-averaged column densities are given in Table [H and 



the corrected source-averaged column densities in Table|5] In the 
remainder of the paper a '^C/'^C ratio value of 53 for the 4 kpc 
molecular ring has been used (Wilson & Rood, 1994). Note that 
this value is uncertain, however, and known to depend on the 
distance to the galactic center However, for simplicity we have 
assumed the same ratio for all sources. 

The rotation diagram results by Ivan der Tak et alj (l2000h 
were not corrected for optical depth, because no isotopo- 
logues were observed. In our data, however, lines for Hj^CO 
and H2C'^0 are observed. The optical depth of H2CO has 
been determined from the 3i2-2ii transition. The observed 
H2CO/Hi^^CO ratio ranges from 6.5 for NGC 6334 IRSl to 37 
for AFGL 2591 and NGC 7538 IRSl, lower than the expected 
ratio, which implies that the main isotopologue is optically thick 
in many of its transitions. The results fo r W 3(H20) are con- 
sistent with the optical depth found by He lmich & van Dishoecki 
(Il997h . The optical depth of H'^CO was estimated based on the 
52,3- 42,2 transition observed for both H^^CO and H2C"*0. The 
H2^CO/H2C'**0 ratio is lower than the expected value of 5-8 for 
most sources, but it is unclear whether this difference is signifi- 
cant. We assume therefore for simplicity that H'^CO is optically 
thin for all sources. 

No direct comparison of the same transition is possible for 
CH3OH with '^CHsOH. Many detected lines have strengths that 
are lower by more than an order of magnitude compared to the 
=5trongest transitions. It is therefore possible to determine the col- 
umn density of the main isotopic species by excluding lines pre- 
;2jdicted to be optically thick for the source-averaged column den- 
.2) sity from the rotation diagram. This also circumvents the addi- 
.2)tional error introduced by the uncertain isotopic ratios. As seen 
■ l)in Fig. lA.2l this excludes mainly lines lying below the fit to the 
■2) optically thin lines and results in a much better determined rota- 
■2)tion temperature. 

The limits on the HNCO/HN'^CO ratios are smaller than 
~53, and HNCO could thus be optically thick for all sources. 
Note that the optical depth determination of HNCO is poten- 
tially complicated by the close spacing of the lines of the main 
species and its isotopologue of only ~5 MHz. However, since 
most of our sources have line widths of only 3-4 km s the 
upper limits for HN'^CO are significant. Since HN'^CO has not 
been detected for any of our sources and most sources studied by 
iZinchenko et alj (I2OOO1) . HNCO is assumed to be optically thin 
throughout the remainder of this paper. 

Since it was not one of the target molecules, CH3CN was 
not observed for all sources. Its lines and those of CH^^CN 
were coincidentally present in the image sidebands for some of 
the sources. For G24.78, NGC 6334 IRSl, and W 33 A, both 
the main species and its isotopologue have been detected and 
CH3CN is found to be optically thick in the 12-11 and 13-12 
transitions with ratios ranging from 5.8 for G24.78 to 24 for 
N GC 6334 IRSl For G24.78, this result confirms the findings 
bv lBeltranetal.l(l2005h . 

Since HCOOCH3 and CH3OCH3 have relatively high col- 
umn densities, predictions were made for which lines were opti- 
cally thick. All lines detected for both molecules were found to 
be optically thin, if we assume that the emission does not come 
from a radius smaller than ^r=i()0K- 



4.3. Rotation temperatures 

In the following subsection we discuss each molecule individ- 
ually and compare it with literature; the main results are sum- 
marized at the end of this subsection. The rotation temperatures 
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Fig. 3. Spectra centered at 240.25 and 353.05 GHz for all sources. Note that the Tmb scales are different for individual sources. 
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Table 4. Column density or line intensity ratios for HjCO/H^^CO, H'^CO/HzCi'^O, HNCO/HN'^CO, and CHaCN/CH'^CN before 
correction for optical depth. 



Species 


AFGL 2591 


G24.78 


G75.78 


Sources 
NGC 6334 IRSl 


NGC 7538 IRSl 


W 3(H20) 


W33A 


HjCO/H^'CO" 


37 


9.3 


22 


6.5 


37 


31 


14 


Hi^CO/HaC'^'O* 


<2.4 


2.8 


>1.7 


2.4 


<1.6 


>3.0 


>3.5 


HNCO/HN"CO' 


>2.4 


>11 


>6.5 


>8.6 


>2.0 


>3.6 


>7.3 


CHjCN/CH^^CN'' 




5.8 


>5.9 


25 


7 


>9.5 


7.6 



" Derived from the 3i,2-2i,i transition. 
* Derived from the 52,3-42,2 transition. 

Derived from the 10o,io-9o,9 transition. 
'' Derived from the overall column density determined from all detected transitions. 

Table 5. Column densities for hot core molecules (source-averaged) and cold molecules (beam-averaged) " [cm^^]. 



Sources 



Species 


AFGL 2591 


G24.78 


G75.78 


NGC 6334 IRSl 


NGC 7538 IRSl 


W 3(H20) 


W33A 


Hot molecules 


H2CO'' 


1.3(16) 


6.4(16) 


2.2(16) 


1.3(17) 


2.5(16) 


1.8(17) 


5.4(16) 




2.4(14) 


1.2(15) 


4.0(14) 


2.5(15) 


4.8(14) 


3.3(15) 


1.0(15) 


H2C"*0 


7.4(13) 


4.3(14) 


<2.5(14) 


1.1(15) 


1.3(14) 


<1.7(15) 


<3.0(14) 


CHjOff 


4.7(16) 


2.8(17) 


1.1(17) 


9.7(17) 


1.2(17) 


1.0(18) 


2.0(17) 


'^CHiOH 


<1.8(15) 


2.9(16) 


4.3(15) 


8.5(16) 


6.7(15) 


3.2(16) 


6.5(15) 


C2H5OH 


<1.0(15) 


7.1(15) 


<2.4(15) 


1.9(16) 


5.7(15) 


8.4(15) 


4.7(15) 


HNCO 


8.3(14) 


5.0(15) 


3.8(15) 


4.3(15) 


2.3(15) 


4.9(15) 


6.6(15) 


HN'^CO 


<4.7(14) 


<9.3(14) 


<1.2(15) 


<3.8(14) 


<1.2(14) 


<2.7(15) 


<1.2(15) 


NH2CHO 


<1.3(15) 


7.2(14) 


2.0(14) 


7.4(14) 


5.7(14) 


1.3(15) 


2.1(15) 


CH3CN 


<3.5(15) 


5.9(16)'' 


1.8(15) 


2.9(16)'' 


<8.2(15) 


7.0(15) 


2.7(16) 


CHf CN 


<5.3(14) 


1.1(15) 


<3.0(14) 


5.6(14) 


<1.2(15) 


<7.2(14) 


<5.0(14) 


C2H5CN 


<7.5(14) 


4.0(15) 


<1.2(15) 


5.1(15) 


<9.2(14) 


4.5(15) 


<2.1(15) 


HCOOCH3 


<2.4(16) 


3.1(16) 


7.1(15) 


1.2(17) 


1.4(16) 


5.2(16) 


2.5(16) 


CH3OCH3 


<7.7(15) 


1.2(17) 


2.3(16) 


5.8(17) 


<1.6(16) 


1.5(17) 


2.7(16) 


Cold molecules 


CH2CO 


1.6(13) 


2.1(14) 


6.4(13) 


7.2(14) 


9.7(13) 


1.1(14) 


6.3(13) 


CH3CHO 


3.1(12) 


<1.5(13) 


<2.1(13) 


<1.2(14) 


2.8(13) 


3.5(13) 


<3.0(13) 


HCOOH 


3.4(13) 


2.1(14) 


7.9(13) 


4.9(14) 


9.8(13) 


1.5(14) 


1.3(14) 


CH3CCH 


7.9(14) 


2.3(15) 


8.6(14) 


5.2(15) 


8.4(14) 


1.5(15) 


1.3(15) 



" The column densities for the hot core molecules are corrected for beam-dilution using the values of z/bf listed in Table|2] The cold molecules are 
given as a beam-averaged column density for a beam of 14". 
* From CO multiplied by 53. 

From rotation diagram excluding optically thick lines. 
'' From CHf CN multiplied by 53. 



were derived from single temperature fits to rotation diagrams 
(see Sect. 13.11 ). Calibration uncertainties of 15% are shown as 
error bars on the rotation diagrams. All rotation diagrams are 
shown in Appendix [A) and resulting rotation temperatures in 
Table |6l 

H2CO: For H2CO, the data are taken from' van der Tak et alJ 
(HoOOrexcept for G24.78 and G75.78 (see Fig. \K^. The ro- 
tation temperatures for all sources except NGC 6334 IRSl and 
W 3(H20) are -80 K. NGC 6334 IRSl and W 3(H20) have tem- 
peratures of about 100 K higher However, as the optical depths 
are similar for most sources, this temperature difference between 
the sources is real. 

CH3OH: For CH3OH no clear temperature trend is seen, but 
the temperature is generally higher than 100 K. This implies that 
the emission is coming from warm gas. The CH3OH rotation 
diagram for some sources like W 33A appear s to consist of two 
temp erature components as seen previously bv lvan der Tak et all 
(I2OOO1) (Fig. |A.2| i. However, deducing the presence of two com- 



ponents from a rotation diagram is non-trivial as this behavior 
could also result from sub-thermal excitation and optical depth 
effects. Indeed, excluding lines of high optical depth leads to the 
disappearance of the two-component structure for most sources 
other than W 33A. The fits in Fig. lA.2l are especially sensitive to 
lines with high excitation temperatures and are therefore thought 
to probe the hot component of the CH3OH gas. 

CH2CO: CH2CO has temperatures that are very low, gener- 
ally only around 40 K. Slightly higher temperatures are found 
for NGC 6334 IRSl and W 3(H20). These temperatures sug- 
gest that this molecule is not present in hot regions, but rather 
in the cold envelope. The IRAM and JCM T detections for 
NGC 7538 IRS 1 and W 3(H20) are consistent. MacDonald et all 
(119961) find very high values of 273 K for G34.3 based on three 
lines with /Su >100 K. However the integrated intensity for the 
CH2CO 17i,i5-16i,i5 lin e at 346.6 GHz wi th >150 K in this 
survey, also detected by iMacDonald et al.l (Il996l) . is consistent 
with the low rotation temperature found in our survey. 
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Sources 








Species 


AFGL2591 


G24.78 


G75.78 


NGC 6334 IRSl 


NGC 75 3 8 IRSl 


W 3(H20) 


W33A 


H2CO 


89" 


83±9 


69+6 


193" 


87" 


181" 


88" 


CH3OH 


147±11 


211±13 


113±7 


178+10 


156+10 


139+8 


259+16 


C2H5OH 


[139] 


104±7 


[139] 


166+15 


164+17 


129+16 


122+16 


HNCO 


64+8 


96±11 


108+13 


105+12 


278+88 


147+24 


85+8 


NH2CHO 


[119] 


170±20 


78+22 


166+29 


164+28 


71+7 


40+13 


CH3CN 


[218] 


261 ±20 


189+25 


170+13 


[218] 


196+14 


278+44 


CzH^CN 


[96] 


99±5 


[96] 


92+3 


[96] 


94+16 


[96] 


HCOOCH3 


[119] 


121±6 


87+13 


144+7 


134+8 


109+7 


112+15 


CH3OCH3 


[130] 


134±12 


[130] 


241+35 


[130] 


94+5 


43+7 


CH2CO 


[47] 


39±6 


40+6 


67+20 


37+4 


54+8 


45+7 


CHiCHO 


[20] 


[20] 


[37.5] 


[37.5] 


18+5 


16+3 


[37.5] 


HCOOH 


[55] 


42±7 


[55] 


75+15 


73+12 


189+108 


38+6 


CH3CCH 


45±12 


59±6 


59+14 


63+12 


58+10 


82+15 


40+6 



" The bracket not ation indicates assumed r otation temperature, when no temperature could be derived from a rotation diagram. 
" H2CO data from lvan der Tak et alj ilOM) . 



CH3CHO: Due to unfavorable weather conditions, IRAM 
observations of CH3CHO have only been performed for 
AFGL 2591, G75.78, NGC 7538 IRSl, and W 3(H20). 
Almost all targeted lines were detected toward the sources 
NGC 7538 IRS 1 and W 3(H20). Both sources have r otation tem- 
peratu r es <20 K, consistent with the results found bv llkeda et al.l 
(I2OOIL I2OO2I) . Although its relatively high dipole moment of 
2.69 D implies that sub-thermal excitation can play a role, it is 
unlikely that the temperature is above 100 K, especially since 
molecules with similarly high dipole moments such as HNCO 
and CH3CN have much higher rotation temperatures. 

C2H5OH: C2H5OH consistently has temperatures of more 
than 100 K. Thi s is somewhat higher than was found by 
llkeda et al.l(l200lb . which could be due to both the detection of 
more high excitation lines in this study. For the lower energy 
transitions, a steeper slope is observed, which indicates a lower 
temperature. 

HCOOH: HCOOH has relatively low temperatures rang- 
ing from -40 K for G24.78, and W 33A to -70 K for 
NGC 6334 IRS 1 and NGC 7538 IRS 1 . For NGC 7538 IRS 1 and 
W 3(H20), it was possible to combine results from both IRAM 
and the JCMT For NGC 7538 IRSl, low rotation temperatures 
are confirmed, whereas for W 3(H20) a scatter of more than 
one order of magnitude is present. As the W 3 (OH) region is 
close to W 3(H20) it is possible that emission from this region 
is picked up in one of the settings. Alternatively, HCOOH could 
be present in both hot and cold gas with only the cold compo- 
nent present for NGC 7538 IRSl. Due to the larger beam size 
of the JCMT, the 230 GHz observations are more sensitive to 
extended emission, whereas the IRAM 270 GHz observations 
probe a more compact region. Since most sources have low rota- 
tion temperatures, we conclude that it is mostly located in a cold 
environment. Low rotation temperatures have previously been 
inferred for HCOOH in other high mass star forming regions 
(llkeda et all l200Tl) . 

HNCO: The rotation temperature of HNCO varies between 
the different sources. This is not unexpected since the HNCO ex- 
citation is known to be do minated by radiation and not collisions 
jChurchwell et aI.L Il986l) . It is therefore not likely to provide 
much information about the actual kinetic temperature of the 
species. However it does need to be close to an infrared source 
to be excited, and thus must be present in hot gas. Finally, it is 
important to note that excluding the line with >300 K (see 



Fig. |A.7t strongly reduces the inferred rotation temperature for 
NGC 6334 IRSl from 180 K to 106 K. 

NH2CHO: NH2CHO has a variety of rotation temperatures. 
The scatter in the derived rotation temperatures appear to be 
relatively large, with especially NGC 7538 IRSl being an out- 
lier (Fig. IA.8b . There is also some indication that there may be 
two temperature components - most clearly seen for the sources 
where many lines with >200 K are detected. Since NH2CHO 
has lines very close to the detection limit, this component may 
also be present for the other sources, but simply below the de- 
tection limit. 

CH3CN: Rotation temperatures of CH3CN range from 170- 
280 K. It is thus always present in hot gas (see Fig. IA.9I ). The 
particularly high temperatures for G 24.78 and W 33A are up- 
per limits due to high optical depths. Pankonin et alj ( 1200 Ih find 
lower rotation temperatures of 89 K for G24.78 which con- 
firm this picture. For NGC 6334 IRSl only three lines with 
£11 > 500 K have been detected, making the resulting rotation 
temperature uncertain. However, it seems clear that the rotation 
temperature for CH3CN is high, and the species must be present 
in hot gas. 

C2H5CN: The other cyanide, C2H5CN, only has enough de- 
tected lines for three of the sources to allow a determination of 
the rotation temperature (see Fig. lA.lOb . For all three sources 
temperatures of ~90 K are found. Fewer lines have been de- 
tected for W 3(H20), but the fact that the results for this source 
are rather similar gives additional credibility to the derived tem- 
peratures. 

CH3CCH: CH3CCH generally has rotation temperatures be- 
low -80 K- except for G24.78, NGC 6334 IRS 1, and W 3(H20). 
The higher temperatures for these three sources are, however, 
questionable. Removing the point with E^ >300 K gives signifi- 
can dy low er temperatures for G24.78 and NGC 6334 IRSl (see 
Fig. lA.llI ). The W 3(H20) generally has a higher uncertainty 
and scatter on the data, which may in fact cause the higher tem- 
perature found for this source. 

HCOOCH3: HCOOCH3 is detected in almost all of the 
sources, but with a relatively large scatter on the data. The 
temperature generally seems to be between 100-150 K, with 
somewhat lower values when less lines have been detected. For 
G75.78 the results are inconclusive as the upper limits on high 
energy transitions are not significant. Somewhat lower tempera- 
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tures were found bv llkeda et alj (1200 Ih . which could be due to a 
difference in the number of detected lines. 

CH3OCH3: CH3OCH3 also appears to have a relatively large 
spread in rotation temperatures (see Fig. IA.13T l. For W 33A, only 
lines with <100 K were detected, which result in low rotation 
temperatures. The upper limits on lines with E^, >150 K are not 
significant, however, and do not exclude temperatures above 100 
K. Rotation temperatures abo ve 100 K are commonly seen for 
CH3OCH3 dIkedaetal.Ll200lh . 

In summary, the temperatures for most molecules are rather 
consistent between all sources. The molecules can clearly be 
divided into "hot" (T > 100 K) and "cold" (T <100 K) 
where HjCO, CH3OH, C2H5OH, HNCO, NH2CHO, CH3CN, 
C2H5CN, HCOOCH3, and CH3OCH3 belong to the first and 
CH2CO, CH3CHO, HCOOH and CH3CCH to the second cat- 
egory. 

4.4. Column densities 

Tableland Figure |4] present the source-averaged column den- 
sities, i.e., corrected for beam dilution (Sect. |3.2| i, for the hot 
molecules per source. Absolute column densities are compared 
to determine the variability between diff'erent sources. Since we 
have assumed the beam dilution to be the same for all molecules, 
the trends for the beam-averaged column densities are similar 
Molecules with rotation temperatures below 100 K are present 
in a different, cold region, presumably the envelope surrounding 
the hot region, and beam-averaged column densities are there- 
fore given for those. Thus, the column densities for the hot and 
cold molecules in Fig. |4]and|5]cannot be compared directly. 

Clearly, most molecules have column densities that vary by 
about an order of magnitude in the hot gas between the different 
sources. The column density variations do not scale with the lu- 
minosity and thus total mass in the hot region. CH3OH, H2CO, 
but also HCOOCH3, and CH3OCH3 - if detected - have col- 
umn densities that are orders of magnitude higher than those of 
C2H5OH, HNCO, and NH2CHO. As will be discussed in Sect.|5] 
abundances show a similar trend. The advantage of comparing 
column densities is that they come directly from the observa- 
tions and no assumption on the H2 column in the hot gas enters 
the analysis. 

Most hot molecules have source-averaged column densi- 
ties that are highest for the line-rich sources NGC 6334 IRSl, 
G24.78, and W 3(H20). For the oxygen-bearing molecules 
C2H5OH, CH3OCH3, HCOOCH3, G24.78 has a higher column 
density than W 3(H20). W 3(H20) has a higher column density 
of CH3OH and H2CO, however Thus H2CO and CH3OH have 
variable column densities compared to the other oxygen-bearing 
molecules. Nonetheless the column density trends between the 
oxygen-bearing species are similar 

The HNCO and NH2CHO column densities peak for W 33A 
and G24.78, and have very similar column density trends with 
respect to each other W 33A has the highest and G75.78 the 
lowest column densities. NGC 7538 IRSl has a relatively high 
column density for these species compared with the oxygen- 
bearing species. The column densities of the nitrogen-bearing 
species do, therefore, behave different from the oxygen-bearing 
species. 

CH3CN has a relatively high column density for G24.78 and 
W 33A, but the uncertainties make it more difficult to compare 
its column density trends with other molecules (see Sect. 14.31 ). 
The general CH3CN column density trends appear more similar 
to those of HNCO and NH2CHO than to the O-bearing species. 
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Fig. 5. Histogram of molecular beam-averaged column densities 
of the cold molecules for all sources. 



However, if we consider only the two sources for which both the 
main species and its isotopologue CH^^CN were detected with 
£11 <500 K, the relative column densities appear more similar to 
those of the oxygen-bearing species. 

Figure |5] shows the beam-averaged column densities for the 
cold molecules, CH2CO, HCOOH, and CH3CCH. Clearly all 
three have very similar trends. NGC 6334 IRSl has the high- 
est column density and AFGL 2591 the lowest. For CH3CHO, 
the number of sources with detections is limited and it was thus 
not possible to compare its column density trends with the other 
species. 

In summary, the molecular abundances vary about an or- 
der of magnitude between the sources. Furthermore, all oxygen- 
bearing species have similar abundance trends, but the abun- 
dances of HNCO and NH2CHO behave differently. 



5. Column density and abundance correlations 

Molecular column densities and abundances can be compared 
between sources in different ways, two of which will be pre- 
sented here. Abundances can be estimated relative to column 
densities of (i) H2 calculated on the basis of CO, (ii) H2 derived 
from dust continuum observations (see Sect. 13.2b or (iii) one of 
the observed molecules, typically CH3OH. The disadvantage of 
method (i) is that CO gas is likely extended over much larger 
scales than the hot molecules due its lower desorption tempera- 
ture. Observations of high excitation lines of CO isotopologues 
which probe only hot gas, will make this method possible in the 
future. In this paper, we use (ii) and (iii), and adopt (i) only for 
the cold molecules. 

Abundances with respect to H2 column densities from mod- 
eling of dust continuum observations (Sect. |T2] l and to CH3OH 
column densities (Sect. 15.2b have the advantage that they allow 
comparisons to "only material at high temperatures" or "species 
probing the regions solely of interest". Both approaches also 
have disadvantages, however Since H2 is not directly observed, 
one has to rely on other tracers, in our case dust. Extrapolations 
of density and temperature profiles are made to scales smaller 
than actually probed by the observations. CH3OH is directly ob- 
served, but this method implicitly assumes that the abundance of 
CH3OH is constant from source to source. If this is not the case 
incorrect correlations of abundances with respect to CH3OH can 
be found, as we will return to in Sect. 15.21 
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Fig. 4. Histogram of molecular column densities of hot molecules corrected for beam dilution for all sources for a source size 
defined by T >100K. 



CH OH 




C H OH CH OCHO 

2 5 3 

Fig. 6. Empirical correlations between the different molecules 
based on the Pearson correlation coefficients. All solid lines in- 
dicate correlations coefficients >0.8 and the thick lines >0.9. 
When two species correlate, it does not mean that they are di- 
rectly chemically related. Two species can also correlate if they 
are both formed from the same mother species. 



5.1. Abundances inT >100 K gas 

Table [T] Usts the molecular abundances within the gas with 
T >100 K for each source in our sample calculated through 
Eq. Q. The abundances have the same trends between the differ- 
ent sources as the column densities. Most molecules have abun- 
dances that vary one order of magnitude from source to source, 
as was also seen for CH3OH by van der Tak et al. (2000). The 
standard deviation on the average abundances is larger than our 
estimated uncertainty of a factor 2-3 (Sect.|3]l. Thus the observed 
variations are real. Furthermore, the overall abundances are rel- 
atively high compared to values in the literature. This is because 
we have determined the abundances with res pect to the H2 col- 
umn density within ^7-=ioo k, whereas in e.g.. llkeda et al.l (1200 lb 
they are determined based on CO transitions, which trace gas 
that is cooler and extended with respect to the hot core. 

The abundances of the species detected in our sources with 
T > 100 K were compared by estimating the Pearson correla- 
tion coefficients between log(X/H2) vs log(Y/H2) (see Table |8]l. 
C2H5CN is not compared with the other species due to the 
low number of sources with detections. Although the number 
of sources for the other molecules are limited and correlations 
therefore only suggestive, a number of interesting trends are 



seen. These result from groups of molecules with very similar 
abundance trends. 

First, the abundances of the nitrogen-bearing species do not 
correlate with those of oxygen-bearing species. Strong correla- 
tions are, however, found between the oxygen-bearing species 
H2CO, CH3OH, C2H5OH, HCOOCH3, and CH3OCH3. FigurelT] 
shows the abundances of all hot oxygen-bearing species with 
respect to H2 compared to CH3OH/H2. The nitrogen-bearing 
species HNCO and NH2CHO also appear to be well corre- 
lated, which follows from the column density trends (Sect.|4]and 
shown in Fig.|8]i. The NH2CHO abundance for G75.78 is uncer- 
tain, as only few lines were detected. The correlation coefficient 
is 0.38 with this source included and 0.92 without. 

If the HCOOH emission is assumed to arise in hot instead 
of cold gas, its abundance correlates with very few molecular 
species. Its strongest correlation is with HNCO with a coeffi- 
cient of 0.98. This result is surprising and due to the low rotation 
temperatures for HCOOH questionable. 

The abundances of the hot molecules presented in Table |7] 
and Fig. [T] assume the same beam dilution factor correspond- 
ing to 7? 100 K- As Table |6] shows, the rotation diagram tempera- 
tures of several species are significantly higher than 100 K. The 
T -200 K radius is typically a factor of ~3 smaller than that for 
100 K so that the beam dilution and abundances can be a factor 
~9 larger If all emission of hot molecules would arise from a 
smaller region, the correlations would remain the same. If some 
species are present in 100 K gas and others in 200 K gas, the 
correlations would disappear, however. The differences in rota- 
tion temperatures in Table |6] are assumed to arise largely from 
other effects such as optical depth, sub-thermal excitation, or the 
detection of very high or low lying energy levels affecting the 
overall fit. Future interferometric observations of multiple high 
mass sources could determine whether the assumption of a sin- 
gle source size for all molecules is correct. 

For the molecules that have abundances that do correlate, 
the abundance ratios were calculated with respect to either 
CH3OH/H2 or HNCO/H2 and are shown in Table |9] Figure |6] 
shows the empirical relations for molecules with Pearson coef- 
ficients of 0.8 and higher Per set of species the ratios vary no 
more than a factor of three, the previously found uncertainty in 
the abundances. H2CO, HCOOCH3, and CH3OCH3 have 5-10 
times lower abundances than CH3OH. NH2CHO is a factor 3 
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Table 7. Molecular abundances iN(X)/N(li2))- 





AFGL2591 


G24.78 


G75.78 


NGC 6334 IRS 1 


NGC 75 3 8 IRS 1 


W 3(H20) 


W33A 


Average 


cr" 


Hot molecules'^ 


H2CO 


1.7(-7) 


1.6(-7) 


1.8(-7) 


5.4(-7) 


I.2(-7) 


I.0(-6) 


2.1(-7) 


3.4(-7) 


3.2 


CH3OH 


7.0(-7) 


6.5(-7) 


9.2(-7) 


4.0(-6) 


5.7(-7) 


5.4(-6) 


7.3(-7) 


1.9(-6) 


2.0 


C2H5OH 


<1.4(-8) 


1.9(-8) 


2.0(-8) 


7.5(-8) 


2.5(-8) 


4.4(-8) 


1.8(-8) 


3.6(-8) 


2.3 


HNCO 


7.0(-9) 


1.3(-8) 


3.2(-9) 


1.8(-8) 


I.I(-8) 


2.7(-8) 


2.5(-8) 


1.6(-8) 


0.91 


NH2CHO 


<1.7(-8) 


1.8(-9) 


1.7(-9) 


3.1(-9) 


2.6(-9) 


5.6(-9) 


8.5(-9) 


3.9(-9) 


2.7 


CH3CN 


<4.6(-8) 


1.5(-7) 


1.5(-8) 


1.2(-7) 


<3.9(-8) 


3.9(-8) 


1.0(-7) 


8.5(-8) 


5.6 


C2H5CN 


<9.9(-9) 


1.0(-8) 


<1.0(-8) 


2.1(-8) 


<4.4(-9) 


2.5(-8) 


<8.I(-9) 


1.7(-8) 


0.77 


HCOOCHi 


<3.2(-7) 


7.3(-8) 


5.9(-8) 


4.6(-7) 


6.7(-8) 


2.9(-7) 


9.6(-8) 


1.7(-7) 


1.6 


CH,OCHi 


<1.0(-7) 


3.0(-7) 


1.9(-7) 


2.4(-6) 


<7.6(-8) 


8.3(-7) 


1.0(-7) 


7.7(-7) 


9.6 


Cold molecules'' 


CH2CO 


3.1(-11) 


8.8(-10) 


8.5(-ll) 


2.8(-10) 


4.1(-10) 


I.3(-10) 


4.1(-I0) 


3.2(-I0) 


2.9 


CH,CHO 


6.1(-12) 


<6.3(-ll) 


<2.8(-ll) 


<4.7(-ll) 


1.2(-10) 


4.1(-11) 


<I.9(-10) 


5.6(-Il) 


5.8 


HCOOH 


6.7(-ll) 


8.8(-10) 


I.O(-IO) 


1.9(-10) 


4.2(-10) 


I.8(-10) 


8.4(-I0) 


3.8(-I0) 


3.4 


CH3CCH 


1.5(-9) 


9.6(-9) 


1.1 (-9) 


2.1 (-9) 


3.6(-9) 


1.8(-9) 


8.4(-9) 


4.0(-9) 


3.5 



" <r refers to the same order of magnitude as the average. 

* The abundances are derived with respect to the column density o f H2 within Rt=ioo k in S ect. |3.2| and Table[3l 

Th e abundances are deriv ed u sing the CO co l umn d ensities from lvan der Talc et alj i2000h , except G24.78 and G75.78 for which the data comes 
from lCesaroni et 

Ml ( l2003h and lHatchell et alj ( Il998h respectively. 
a refers to the same order of magnitude as the average. 



Table 8. Abundance correlations with respect to H2 for the abundances within Rt=iqo k- for hot core molecules. 





H2CO 


CH3OH 


C2H5OH 


HNCO 


NH2CHO 


CH3CN 


HCOOCH3 


CH3OCH3 


HCOOH 


H2C0 


1.00 


0.98 


0.76 


0.46 


0.53 


-0.07 


0.88 


0.75 


0.76 


CH30H 


0.98 


1.00 


0.86 


0.38 


0.39 


-0.07 


0.92 


0.86 


0.70 


C2HSOH 


0.76 


0.86 


1.00 


0.02 


0.22 


0.12 


0.92 


0.95 


0.30 


HNCO 


0.46 


0.38 


0.02 


1.00 


0.38 


0.81 


0.16 


-0.27 


0.98 


NH2CHO 


0.53 


0.39 


0.22 


0.38 


1.00 


0.17 


0.42 


-0.02 


0.65 


CH3CN 


-0.07 


-0.07 


0.12 


-0.81 


0.17 


1.00 


0.29 


0.20 


-0.51 


HCOOCH3 


0.88 


0.92 


0.92 


0.16 


0.42 


0.29 


1.00 


0.90 


0.49 


CH3OCH3 


0.75 


0.86 


0.95 


-0.27 


-0.02 


0.20 


0.90 


1.00 


0.09 


HCOOH" 


0.76 


0.70 


0.30 


0.98 


0.65 


-0.51 


0.49 


0.09 


1.00 


" Assuming that the HCOOH emission comes from within Rt=ioo k- 



Table 9. Average abundance ratio for our sources and other well-known high mass YSOs with respect to CH3OH for the oxygen- 
bearing species and with respect to HNCO for NH2CHO. cr is the standard deviation between the sources in our sample. 



Molecule 


abundance ratio(cr) Orion Compact Ridge" 


G34.3'' 


Sgr B2(N)'- 


0327.3" 


GC clouds'- 


with respect to CH3OH 


Ho CO 


0.22±0.05 


0.10 


0.009 


>0.0025 




-0.01 


CoH^OH 


0.025±0.013 


0.004 


0.19 


0.005 


0.05 


-0.04 


HCOOCH3 


0.10±0.03 


0.08 


0.86 


0.005 


0.96 


-0.04 


CH3OCH3 


0.31+0.20 


0.05 




0.015 


7 


-0.04 


with respect to HNCO 


NH2CHO 


0.36±0.08 






0.23 


0.07 





"'Sutton et al 

* MacDonald et al 



lll995h 



■199^) 

Nummelin et al. ( 2000h. 



'' pTbbeta l. (2000a). 
''iReauena-Torres etaP j2006h. 



times lower than HNCO. Possibly, coefficients between 0.8-0.9 
are present between species that are not directly chemically re- 
lated, but that are both strongly correlated with a third species. 



5.2. Abundances with respect to CH^OH 

The Pearson correlation coefficients of log(X/CH30H) vs. 
log(Y/CH30H) differ significantly from the correlations in 
Sect. 15.11 First, strong correlations are not present between the 
same species. This is due to the variability of the CH3OH abun- 
dance between the different sources, and to the dependence of 
the correlations on three (X, Y, and CH3OH) rather than two 
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[CH^OH/HJ 



Fig. 7. Comparison of molecular abundances with respect to H2 within /?r=ioo k for the oxygen-bearing molecules. The gray lines 
indicate the linear fit to the data. An error-bar representing the calibration uncertainties is depicted in the lower right corner 




[HNCO/H ] 



Fig. 8. Comparison of molecular abundances with respect to H2 
for HNCO with NH2CHO within Rt=wok- The source indicated 
in gray is G75.78 and the limit is AFGL 2591. An error-bar rep- 
resenting the calibration uncertainties is depicted in the lower 
right corner. 



(X and Y) species. The abundances of two species with respect 
to CH3OH will correlate when their abundances with respect to 
H2 correlate but also when the abundance of one of the species 
correlates with that of CH3OH. This makes it difficult to use cor- 
relations of abundance ratios with respect to CH3OH to chem- 
ically link molecules. Instead it is better to compare X/CH3OH 



with the CH3OH column density or abun dance, i.e. N(CH30H) 
or CH 3OH/H2 as has been done by e.g., iRequena- Torres et akl 
(120061) . 



5.3. Comparison to other well known sources 

The abundance ratios with respect to CH3OH are compared in 
Table |9] to other well-studied high mass star forming regions. 



specifica lly the Orion C ompact Ridge (Sutt on et al.l 1995h. 



G34.3 (iMacDonald et al.L ir99&i. Ser B2(N) (iNummeUn et all 
l2000h . and G327.3 (iGibb et al.Ll20"00ah . The ratios for the Orion 



Compact Ridge all agree within the error bars for H2CO and 
HCOOCH3. CH3OCH3 has a relatively low abundance with re- 
spect to CH3OH. Only three lines have been detected for this 
species in the Compact Ridge and it has a rather low excita- 
tion temperature of 75 K implying that it is perhaps present in a 
colder region. 

In contrast, the relative abundances for G34.3 deviate 
more from those in ou r work. The H2CO/CH3OH ratio by 
iMacDonald et al.l d 1 9961 ) is probably a lower limit as only one 
line was detected for G34.3 and the column density of H^^CO 
was calculated to be higher than the main H2CO isotopo- 
logue based on their respective line intensities. Both C2H5OH 
and HCOOCH3 have relatively high abundances with respect 
to CH3OH. This could potentially be due to high optical 
depth in the lines of CH3OH. MacDonald et al. ( 1996) assumed 
the emission to be optically thin, and thus underestimate the 
CH3OH column density. The relative abundances of C2H5OH 
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and HCOOCH3 are similar between G34.3 and the sources stud- 
ied in this paper. 

T he abundance rat i os with respect to CH3OH for Sgr B2(N) 
from iNummelin et alj ( l200tf ) appear to be significantly lower 
than for the sources in our sample. This is due to the CH3OH 
emission being decomposed into a core and a halo compo- 
nent for Sgr B2(N), whereas that of the other molecules is not. 
For this source it is more relevant to compare the relative ra- 
tios between the different molecules. The C2H5OH abundance 
is high compared to HCOOCH3 and CH3OCH3. However, the 
HCOOCH3/CH3OCH3 and the NH2CHO/HNCO ratios of 0.33 
are similar to those for our sources^ 

G327.3 in the study of iGibb et al.1 (l2000al) has a 
C2H5OH/CH3OH abundance ratio a factor of 2 higher compared 
to the average for our sample. Small deviations could be caused 
by their relatively low rotation temperature for C2H5OH of 66 K. 
HCOOCH3 and CH3OCH3 have relatively high abundances with 
respect to CH3OH, but HCOOCH3/CH3OCH3 is within a factor 
two of the ratio derived for our sample. Again NH2CHO/HNCO 
is the same within error to our sample. 

Molecular abundance ratios can also be compared to 
those found in shocked gas in the galactic center by 
iRequena- Torres et al. I (l2006h . Their sources have relatively con- 
stant ratios of complex molecules with respect to CH3OH of 
0.04. The H2CO/CH3OH and CH3OCH3/CH3OH ratios are sig- 
nificantly high er for our hot cores. Small di fferences may be due 
to the fact that lRequena-Torres et al.l (12006 ) observed rotational 
transitions with iiu <40-50 K, whereas lines with excitation 
energies that range from 20-900 K have been observed in this 
study. 

In summary, the overall emerging picture is that of relatively 
constant abundance ratios between various molecules compara- 
ble to the sources for our sample. Distinct differences for indi- 
vidual species can often be ascribed to specific observation or 
analysis problems. 

5.4. Comparison with ice abundances 

It is interesting to relate the gas-phase abundances derived 
above to the ice abundances from mid-infrared observations. 
Of the sources in our sample, AFGL 2591, NGC 7538 IRSl 
and W 33 A were also obse rved at mid-infrared wavelengths 
with ISO (Gibb et al.L |2004 . Table [10] lists the abundances of 
the ice species from that paper assuming an H2O abundance 
of 10 with respect to H2. W 33A is particularly rich in or- 
ganics in the solid state whereas the other two sources mainly 
have upper limits. The solid state abundances of CH3OH are 
at least an order of magnitude higher than the comparable 
source-averaged gas phase abundances. The H2CO abundance 
of ~2xl0"^, alth ough more uncertai n, is also larger by an or- 
der of magnitude jKeane et al.Ll200l1) . For W 33 A, however, the 
solid H2CO/CH3OH abundance ratio is consistent with the gas 
phase abundance ratio. 

Where detected toward high mass YSOs, solid HCOOH 
and CH3C HO have abundances of 10~^-10~^ and -10"^ re- 
spectively ('Schutteetafl 119991: lKeaneetal.L I2OOII: iGibb et al.L 
[2000b, 2004). Similar abun dances for s o lid HC OOH are found 
toward background stars bv iKnez et al.l (l2005h . These absolute 
solid state abundances are four to five orders of magnitude higher 
than the gas phase abundances for the same sources. 

The XCN band is detecte d in many obje cts and 
usually assigned to OCN" d Grim & Greenbera, 11987 



Demvketal., 1998; Pendleton et all 
2000b; .van Broekhuizen et al.. ,2005lP oCN~ 



1999; Gibb et al 



is thought to 



convert to and evaporate off grains as HNCO (Sect. 16. It . 
Abundances of solid OCN_vary and range from 1x10"^- 
4x10"^ (van Broekhuiz en et all 2005; Thietal., 2006). Gas 
phase abundances of HNCO are factors of 100-1000 times 
lower. 

Many of the species are ex pected to have abund ances in the 
ice below the detection limit. iBoudin et al. I (Il998h determined 
upper limits for the column densities of a number of more com- 
plex organic molecules, including C2H6 and C2H5OH toward 
NGC 7538 IRS9. The upper Umit on the C2H5OH abundance is 
1.2x10"^ with C2H5OH/CH3OH <0.28. This is consistent with 
the C2H5OH/CH3OH abundance ratio of 0.009 in the gas phase 
for our sample of sources. 

In summary, the absolute ice abundances of observed "first 
generation" species, CH3OH and OCN , are 10-10^ times 
higher compared to the gas phase values, and as much as a factor 
of lO'^-lO^ for HCOOH and CH3CHO. 

6. Discussion and comparison to astrocliemical 
models 

6.1. Generai modei considerations 

Pure gas phase chemistry models such as presented in e.g., 
iLee et al.] (11996'). are unable to explain the high abundances 
of many complex organics, especially of highly saturated 
molecules in hot gas phase environments. CH3OH, for exam- 
ple, is estimated to have abundances of ~10 through pure gas 
phase chemistry, whereas abundances of 10" ^ are found. For 
CH3CN abundances of -10 '° are modeled (" Lee et al.L Il996h 
whereas abundances of ~10"^ are detected. Furthermore, they 
are not able to explain the presence of compl ex organics in ices 
by simple freeze-out of gas phase species (see'Herbst', '2005', for 
a review). For highly saturated molecules, solid state chemistry 
provides a likely alternative. Here one can distinguish basic grain 
surface chemistry and chemistry induced by energetic process- 
ing (UV irradiation, energetic particle bombardment and/or ther- 
mal heating) which can proceed also inside the ice. 

Surface chemistry models have been developed by several 
groups (e.g., iHasegawa & Herbstl 1 1993HCharnlevi 120011: iKeand 
2001) (see Fig. 1). The results are generally very sensitive to 
the atomic C, O, and N abundances in the gas. Moreover, the 
dust temperature (Tdust) plays a crucial role in the mobility of 
the atoms and radicals as well as the availability of CO as a 
reaction partner on the grain, which evapor ates at Tdust >20 K 
(iTielens & Hagenlll982t lCaseni et al.Lll993 l). However, if CO is 
located in an ice environment that is dominated by H2O it will be 
available for reactions in th e ice up to m uch high er tem peratures 
(.Sandford & Allamandoialll993l;ICoUings et al.i.l2003h . 

If grain surface chemistry is the explanation for our observed 
correlations, the constant abundance ratios for different sources 
must imply very similar conditions during ice formation, at least 
of the ice layer that contains the "hot" molecules. This is also 
confirmed by the similarity between the composition of ices ob- 
served toward background stars com pared to those to ward low 
and high mass star forming regions ( Knez et al.L 12005 ). Perhaps 
the simplest explanation is that the bulk of the ice forms at simi- 
larly low temperature conditions in which the bulk of the oxygen 
and nitrogen are in atomic form and C/CO is at a f airly constant 
ratio, around 0.01 as computed for dense clouds dGredel et al.L 
11989 ^. When these "first generation" species evaporate together 
and have similar destruction rates, a constant abundance ratio 
would result even in regions where the temperatures of the gas 
and UV-flux are very different. 
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Species 


W33A'= 


NGC 7538 IRS 1" 


AFGL 2591* 


HH 46'' (low mass) 


CK T (field star) 


gas phase" 




10-5 


10-5 


10-5 


10-5 


10-5 




H2CO 


0.65 










3.2(-7) 


CH3OH 


1.4-1.7 


<0.4 


1.4 


1.2 


<0.11 


2.0(-6) 


HCOOH 


0.37-0.71 


2.2 


<1.7 


<1.5 


0.11 


3.8(-10) 


CH,CHO 


0.98 


1.2 


<0.66 






5.6(-ll) 


OCN- 


0.63 


<0.05 




<0.12 




1.6(-8)^ 



" FromlGibb et alj il2000bh: lDM"tois et al.l i ll999h : ISchutte et alj (ll999h : IWhittet et al.l (l2001h : lDemvk et al.1 ( Il998h . 

' FromGibb et al. (2004). 

From .Knez et al.. ( ,2005) and references therein. 

Average abundances from this paper. 
' Assuming that OCN- evaporates as HNCO. 



In contrast, "second generation" species are expected to peak 
at different times dependent on their formation and destruc- 
tion processes. While this argument has been used to explain 
the anti-coincidences between the nitrogen and oxyg en-bearing 
species ( Caselli et al.Lll993l:lRodgers & Charnleyil2003i) . it is in- 
consistent with our observed correlations of supposedly "sec- 
ond generation" species HCOOCH3 and CH3OCH3 with other 
"first generation" oxygen-bear ing molecules. Indeed, in mod- 
els by iGarrod & HerbstI ( l2006h where HCOOCH3 is formed in 
the solid state, many of the oxygen-bearing species have similar 
gas phase abundance variations throughout the chemical evolu- 
tion and have constant abundance ratios. Solid state formation 
is therefore more likely for the supposedly "second generation" 
oxygen-bearing species. However, the abundances of oxygen- 
bearing and nitrogen-bearing species do not correlate. This can 
be due to a different time-dependence of formation and destruc- 
tion or a variable N/C abundance ratio. 

Laboratory experiments show that beyond a certain mini- 
mum dose of hi'adi ation the reactio n pr oducts often do reach 
a con st ant plateau jGerakines et al], [199 6; Moore & Hudson, 



I2OOOI; Ivan Broekhuizen et all 12004 iHagen et al.. ,197 9) 



Oxygen-bearing species, such as H2CO, and nitrogen-bearing 
species, e.g ., OCN- and NH2CHO, can be formed this way 
dGrim et al.L[l989 ). However, there is not enough information in 
the Uterature on the formation of the more complex molecules 
observed here to make a quantitative test of this scheme. It 
is plausible, though, that an equilibrium between formation 
and destruction is reached in interstellar ices and could be 
responsible for the constant abundance ratios between the 
oxygen-bearing species or HNCO/NH2CHO. 

Recent work by ' Garrod & HerbstI (|2006') combines both 
grain surface chemistry and UV-induced chemistry with thermal 
evolution of the ices during the protostellar phase. In particular, 
radicals created in the ices by ultraviolet radiation become mo- 
bile at the higher dust temperatures during the warm-up phase 
and lead to the formation of many complex organic species. 
This suggests that the ice composition changes with temperature. 
Potentially only the ices on grains that are close to the evapora- 
tion temperature have similar relative abundances to those ob- 
served in the gas phase. Since the detected ice features are ob- 
served in a column, they will contain both cold and warmer ices, 
which complicates the comparison of the ice abundances with 
the gas phase, even within one source. 

Species such as HCOOH and CH3CHO are detected with 
much lower column densities in the gas phase than in the ice 
(see Section l54b and have low rotation temperatures. Possibly, 
these species are formed in the solid state at low temperatures 



and destroyed due to reactions at high tempera tures. Since their 
desorption ternperatures are typically ~100 K dViti et all l2004t 
iBisschop et alilin prepj) . thermal desorption cannot account for 
their low temperatures. However, n on-thermal mechan isms in- 
duced by, for example, cosmic rays ('Leger et"an,ll98 5'). can ex- 
plain the low gas phase abundances if they cause a small fraction 
(10-"*) of the ice to evaporate (see Sect. 16.21 ). Alternatively, cold 
gas phase ion-molecule reactions are a ble to reproduce the abun- 
dances of e.g., CH3CHO and CH2CO dLee etaUlT996l see also 
Sect.lOl- 

Ices are expected to be layered because different molecules 
such as H2O, CO, and N2 form or condense at different den- 
sities. Since nitrogen is transformed into N2 deepest in the 
cloud it is expected to freeze-out last. To form species such as 
HNCO or NH2CHO a more reactive form of nitrogen, i.e. N- 
atoms are needed. N-atoms could potentially be present in a 
non-hydrogen bonding layer dominated by CO or a layer dom- 
inated by H2O. If they were present in the non-hydrogen bond- 
ing layer a small fraction of non-thermal desorption would result 
in cold gas phase N-bearing species. The lack of these cold N- 
bearing molecules and high rotation temperatures of HNCO and 
NH2CHO suggests that they desorb at higher temperatures and 
must thus be present in a different, more tightly bound ice layer 

In summary, pure gas-phase reactions cannot reproduce the 
abundances of most species depicted in Fig.[T] and second gen- 
eration gas phase chemistry is inconsistent with the observed 
correlations for the oxygen-bearing species. Constant abundance 
ratios after ice evaporation can be explained either by similar 
initial conditions for ice formation or an equilibrium between 
formation and destruction due to energetic processes. The low 
abundances of cold gas phase species could be due to gas phase 
ion-molecule reactions or non-thermal desorption of a small 
fraction of the ices e.g., by cosmic rays. 

6.2. Individual molecules 

This subsection elaborates on the general conclusions stated in 
Sect. l6.1l for individual molecules, grouped according to the dif- 
ferent branches of Fig. [T] 

H2CO and CH3OH: CH3OH and H2CO are molecules com- 
monly expected to originate from grain-surface chemistry. They 
are formed through successive hydrogenation of CO on grains, 
a scheme that has been tested experimentally by iHiraoka et al.l 
( l2002h .lWatanabe et al.' ("2004"), and Fuchs et al. (in prep.). Both 
H2CO and CH3OH have been detected in these experiments with 
different abundance ratios dependent on the exact conditions. 
H2CO can also be formed by UV-irradiation of CH3OH con- 
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taining ices (lAllamandola et all Il988t iBernstein etall Il995h . 
However, the high ice abundances of CH3OH (e.g.. lGibb et al.L 
[2004) and the constant abundance ratio between H2CO/CH3OH 
in this paper are inconsistent with this CH3OH destruction mech- 
anism. Our data imply that hydrogenation reactions are more 
feasible formation mechanisms. 

Furthermore as discussed in Sect. 15.41 the solid state abun- 
dance of CH3OH appears to be at least an order of magnitude 
higher than the comparable source-averaged gas phase abun- 
dances. This could be due to a number of effects: i) CH3OH is 
partly destroyed upon evapor ation, ii) intrinsic CH3OH ice vari- 
atio ns from sourc e to sou rce ( Brooke et al. '. 1999; Dartoi s et all 
11999; Pontoppi dan et al.L |2004|) . iii) an underestimate of the 
beam-filUng factor of the CH3OH emitting gas. Possibly, all 
three factors contribute to the relatively low gas phase abun- 
dance. 

CH2CO, CH3CHO, and C2H5OH: In Figure [T] CH2CO, 
CH3CHO, and C2H5OH are thought to be interrelated through 
successive hydrogenation as well as to H2CO and CH3OH by 
successive H- and C-atom addition. All three speci es are ex- 
pecte d to evaporate from the ice together with H2O (IViti et al.L 
I2OO4I) and have similar rotation temperatures. This is not the 
case, however. Indeed, interferometric observations of CH3CHO 
in high mass star forming regions show that CH3CHO emis- 
sion is extended with respec t to hot cores species such as 
HCOOCH3 and CH3OCH3 (S [2005 ). This suggest that 
CH2CO, CH3CHO, and C2H5OH gas are not present in the same 
region. Possibly, H-atom addition to CH2CO and CH3CHO on 
grains is very efficient and completely converts them to C2H5OH 
at high temperatures, so that only C2H5OH evaporates from the 
ice and is detected in the warm inner regions. Non-thermal des- 
orption of a small fraction due to e.g., cosmic ray spot heating of 
CH2CO and CH3CHO ice in cold regions could explain the low 
gas phase abundances and rotation temperatures, as well as the 
high ice abundances of CH3CHO. Alternatively, C2H5OH could 
be formed through another surface reaction related to CH3OH 
and H2CO formation, explaining its constant abundance ratio. 

For CH2CO and CH3CHO there are also potential gas 
phase formation mechanisms. Millar & Nejad ( 1985) argue that 
CH2CO could result from radiative associative reactions be- 
tween CH3 and CO. Assuming that CH3CO^ dissociative re- 
combination results in CH2CO for 50% of the cases, this would 
allow for CH2CO abundances of ~10"^, comparable to our ob- 
served abundanc es and also to tho se by Turner et al. ( 1999). Gas 
phase models bv lLee et al.l (1 19961) predict fractional abundances 
of 5x10 fo r CH3CHO a nd its iso mers. This is lowe r than the 
abundances of llkeda et al] ([2001) andl Charnle"vl (l2004 but close 
to th e abundances for some of our sources. He rbst & LeungI 
([1989) proposed that H3O+ and C2H2 react to form CH3CHO 
after dissociative recombination, which results in higher abun- 
dances of up to ~10"'" in dark clouds. Many of these conclu- 
sions depend on high branching ratios toward these products in 
dissociative recombination. If three body break-up prevails as 
found f or some species, pure gas phase production becomes un- 
tenable ( lGeppertetaUl2005l) . 

HCOOH: High abundances of HCOOH in the ice, low gas 
phase rotation temperatures, and orders of magnitude lower 
abundances in the gas phase compared to the solid state can 
be due to two reasons: i) HCOOH is quickly destroyed in the 
gas phase upon evaporation, and ii) HCOOH is abundant in low 
temperature ice, but destroyed at higher temperatures due to re- 
actions in the ice. Gas phase destruction is a possibility, but it 
is unclear why the destruction rate of HCOOH would be or- 
ders of magnitude larger compared to other species. Even if this 



would be the case, most high mass star-forming regions are ex- 
pected to have an evaporation front that moves outward through 
the envelope as the luminosity increases. HCOOH should thus 
continuously freshly evaporate from icy grains. The second pos- 
sibility of HCOOH destruction in the ice at high temperatures 
is much more likely and can be tested in the laboratory. The 
lower gas phase abundance can then be explained by a small 
fraction of the ice (~10""^) that evaporates due to non-thermal 
desorption mechanisms. The low rotation temperatures of ~40- 
70 K for most sources in our sample are also consistent with this 
picture. Finally, it can explain the high HCOOH abu ndances in 
shocked regions by e.g. Requena -Torres et al.l ((2006), where the 
ice is warmed up and evaporated over much shorter timescales. 

The precise formation of HCOOH through surface chemistry 
still remains a puzzle, however In Figure[T]the formation mecha- 
nisms in the solid state is either through CO or CO2. Laboratory 
experiments indicate that it only evapora tes at very hi g h tem - 
peratures (~ 100 K) similar to H2O tColUngs et all |2004 . 
This is even the case fo r HCOOH di luted in CO or CO2 ices 
(Bisschop et al. in prep.). [Keanel (1200 l l) predicts the formation of 
HCOOH through successive atom addition to CO and finds that 
the HCOOH i ce concentration is propo rtional to the O/O2 ratio. 
Alternatively, lAllen & Robinsoril (|1977() proposed the reaction 
H CO + OH in the solid st ate leading to HCOOH. New models 
by Garrod & H erbst (2006) have HCOOH gas phase abundances 
that are much higher than found for our sources with most of the 
HCOOH actually formed in the gas. If this occurs at intermedi- 
ate temperatures, a fraction of the HCOOH may freeze out on 
grains. Further laboratory experiments are needed to elucidate 
the exact formation mechanism for HCOOH. 

HNCO and NH2CHO: HNCO formation has been explained 
by surface chemistry formation and g as phase reactions 
(Turner et al., 1999t [Keanei l200lh . Both [Turner et al.1 d 19991) 
and Keang (i200ll) argue that NH2CHO is formed on grain- 
surfaces due to its high level of hydrogen saturation. The 
strong correlation between the two species in Sect. [5T| suggests 
they are related. Alternatively, both OCN" and NH2CHO can 
form in ices due to ene r getic processes such as UV-irradiation 
("van Broekhu izen et al.L 



2004t iMufloz Caro & Schuttel 120031: 



Allamandola et al.L 119991) . Van Broekhuizen et al.(2004) con- 
clude that thermal formation of OCN" can, however, also not 
be ruled out. Concluding, solid state processes appear to be the 
most likely formation mechanism for both species, whether they 
involve energetic processing or not. 

HCOOCH3 and CH3OCH3: HCOOCH3 and CH3OCH3 
were previousl y thought to form fr om protona t ed m ethanol in 
the gas phase iMillar et all Il99lh . iHorn et all (1200 4) recently 
found the barrier for the gas phase formation of HCOOCH3 
to be higher than previously anticipated. Moreover, the gas 
phase routes assume that H2COOCH3 and (CH3)20H^ recom- 
bine to HCOOCH3 and CH3OCH3 respectively. Measurements 
by Gepp ert et al.l (l2004t l2005 h have shown that dissociative re- 
combination of CH3OH2 mostly results in a three-body break- 
up instead of its deprotonated counterpart. If this would also 
be the case for species such as HCOOCH3 and CH3OCH3, 
they are more likely to be at least partially formed on grain 
surfaces. The constant abundance ratios between both species 
and CH 3OH also confirm this scenario. Indeed Garrod & HerbstI 
(120061) explain surface formation of HCOOCH3 and CH3OCH3 
through reactions of CH3 O and HCO radicals, first proposed by 
I Allen & Robinsonid 19771) . 

CH3CN, C2H5CN, and CH3CCH: CH3CN can be formed 
from the fast association reaction CHt + HCN and the radia- 
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tive association reaction of CH3 + CN dCharnlev et al.L Il992h . 
The abundances in the gas phase for Orion are, however, ~10 '° 
jCha rnlev et al., 1992) compared to 10 in our sources. The rel- 
atively high abundances plus their high degree of hydrogenation 
for CaHgCN and CH3CN imply that they are possibly formed 
on grain-surfaces. CH3CCH usually has very low rotation tem- 
peratures, and its formation is explained in dark and translucent 
clouds through gas phase ion-molecule or neutral-neutral reac- 
tions. 

In summai-y, most species, such as H2CO, CH3OH, 
C2H5OH, HCOOCH3, and CH3OCH3 probably result from 
"first generation" chemistry in ices, likely dominated by surface 
chemistry. This is possibly also true for CH3CN and C2H5CN. 
CH2CO, CH3CHO, HCOOH and CH3CCH have all been de- 
tected in the gas with low rotation temperatures and, for some 
species, abundances that are orders of magnitude lower than 
those in ices. Non-thermal desorption of a small fraction of the 
ice (~10 is sufficient to explain the abundances in the cold 
gas phase. For CH2CO, CH3CHO, and HCOOH formation plus 
destruction in the ice at higher temperature are plausible, but gas 
phase chemistry cannot be excluded. Models of gas phase re- 
actions are able to accurately reproduce the observed CH3CCH 
abundances. 

7. Summary and Conclusions 

We have performed a survey of selected frequency settings for 
7 high mass YSOs. The survey was targeted at detecting species 
thought to be produced by the grain surface chemistry depicted 
in Fig. [U Other species such as CH3CN, C2H5CN, CH3CCH, 
HCOOCH3 and CH3OCH3 have also been detected. To compare 
abundances between sources, the source luminosity and 850 yum 
dust emission was used to constrain the radius within which T > 
100 K and the column density out to this radius, A^r>io()K- The 
main conclusion derived in this work are: 

- Most detected molecules have similar rotation temperatures 
for all sources. The following molecules can be classified 
as cold (T <100 K): CHjCO, CH3CHO, HCOOH, and 
CH3CCH; and the following as hot (T >100 K): H2CO, 
CH3OH, C2H5OH, HNCO, NH2CHO, CH3CN, C2H5CN, 
HCOOCH3, and CH3OCH3. 

- Groups of molecules have their highest column densities 
for the same sources. This is the case for the oxygen- 
bearing species H2CO, CH3OH, C2H5OH, HCOOCH3 and 
CH3OCH3. HNCO and NH2CHO have their highest column 
in other sources than the oxygen-bearing species, but both 
are highest for the same sources. Due to optical depth effects 
and a limited number of detections it is difficult to establish 
the trends for CH3CN and C2H5CN. 

- Absolute abundances in the hot regions vary by an order 
of magnitude. Relative abundances of the oxygen-bearing 
species and the two nitrogen-bearing species are, however, 
very constant. The absolute abundances are much higher 
than can be produced by gas phase chemistry. The con- 
stant abundance ratios plus the similar rotation tempera- 
tures imply that the oxygen-bearing species, H2CO, CH3OH, 
C2H5OH, HCOOCH3 and CH3OCH3, as well as the two 
nitrogen-bearing species, HNCO and NH2CHO, share a 
common solid state formation scheme. This can be the result 
of very similar conditions during ice formation or a balance 
between formation and destruction in the ice. The oxygen- 
bearing species and the two nitrogen-bearing species are not 
related, however The abundances of two cyanides, CH3CN 



and C2H5CN, are likely due to "first generation" ice chem- 
istry as well. 

- CH2CO, CH3CHO, CH3CCH and HCOOH have very low 
rotation temperatures. The formation of CH2CO, CH3CHO, 
and HCOOH can be explained by solid state or gas phase 
reactions. For HCOOH and CH3CHO gas phase abundances 
are a factor of lO'* lower than those observed in the solid 
state. Solid state formation for these species is thus more 
likely. A low fraction of non-thermal desorption could be 
responsible for the low gas phase abundances. Models of 
ion-molecule reactions in the cold gas phase can, however, 
also reproduce the abundances of CH2CO and CH3CHO. 
CH3CCH on the other hand is probably created by cold gas 
phase chemistry, as its abundance is similar to that found in 
dark clouds. Interestingly, no cold N-bearing molecules are 
found. 

The abundance correlations and branching ratios determined 
in this paper are key to understanding chemical relationships be- 
tween these complex organic molecules and addressing their ori- 
gins. Future spatially resolved observations with interferometers 
such as ALMA will make it possible to confirm molecular re- 
lations and address the chemical variations within each source. 
Additional laboratory experiments can explore the pathways to 
the formation of the most abundant species given these observa- 
tional constraints of constant abundance ratios. Taken together 
such studies will be important to fully understand the chemistry 
of complex organic molecules in star forming regions. 
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Appendix A: Rotation diagrams per molecule for all 
sources 
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Fig. A.2. Rotation diagrams for CH3OH for respectively AFGL2591, G24.78, G75.78, NGC 6334IRS1, NGC 7538IRS1, W3(H20), 
and W33A. The filled black squares are optically thin lines, the open gray squares optically thick lines, based on the estimated 
column densities. The rotation temperatures are fits to the optically thin lines only. 



22 



S. E. Bisschop et al.: Testing grain-surface chemistry in massive hot-core regions 





E (K) 

Fig. A.4. Rotation diagrams for CH3CHO for respectively NGC 7538 IRSl, W3(H20), and W33A. 
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Fig. A.5. Rotation diagrams for C2H5OH for respectively G24, NGC 6334IRS1, NGC 7538IRS1, W3(H20), and W33A. The gray 
triangles indicate upper limits. 
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Fig. A.6. Rotation diagrams for HCOOH for respectively G24.78, NGC 6334IRS1, NGC 7538IRS1, W3(H20), and W33A. The 
gray triangles indicate upper limits. 
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Fig. A.7. Rotation diagrams for HNCO for respectively G24.78, G75.78, NGC 6334IRS1, NGC 7538IRS1, W3(H20), and W33A. 
The filled black squares are the points included and the gray squares excluded from the fit. The gray triangles indicate upper limits. 
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Fig.A.8. Rotation diagrams for NH2CHO for respectively G24.78, G75.78, NGC 6334IRS1, NGC 7538IRS1, WSCHiO), and 
W33A. The gray triangles indicate upper limits. 
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Fig. A.9. Rotation diagrams for CH3CN for respectively G24.78, G75.78, NGC 6334IRS1, W3(H20), and W33A. The gray triangles 
indicate upper limits. 
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Fig. A.IO. Rotation diagrams for C2H5CN for respectively G24.78, and NGC 6334IRS1. The gray triangles indicate upper Utnits. 
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Fig.A.ll. Rotation diagrams for CH3CCH for respectively AFGL2591, G24.78, G75.78, NGC 6334IRS1, NGC 7538IRS1, 
W3(H20), and W33A. The filled black squares are the points included, the gray squares excluded from the fit and the gray tri- 
angles indicate upper limits. 
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Fig. A.12. Rotation diagrams for HCOOCH3 for respectively G24.78, G75.78, NGC 6334IRS1, NGC 7538IRS1, WSCHjO), and 
W33A. The filled black squares are the points included and the gray triangles indicate upper htnits. 
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Fig.A.13. Rotation diagrams for CH3OCH3 for respectively G24.78, NGC 6334IRS1, WSCHiO), and W33A. The filled black 
squares are the points included and the gray triangles indicate upper Umits. 
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Appendix B: Detected lines per species for all 
sources 
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Table B.l. Observed line fluxes J TMsdV (K km s"^) for H2CO and its isotopic species. 



Frequency 


Transition 








Sources 








[GHz] 




AFGL2591 


G24.78 


G75.78 NGC6334IRS1 


NGC 7538 IRS 1 


W 3(H20) 


W33A 


H2CO 


218.222 


3o,3- 2o,2 




17 


14 










218.476 


32,2- 22,1 




5.6 


3.3 










225.698 


3i,2- 2i,i 


15 


27 


22 


57 


31 


34 


22 


363.946 


52,4- 42,3 




15 


5.9 










364.103 


54,1/2- 44,0/1 




7.3 


2.7 










364.275 


53.3-43.2 




16 


7.5 










364.289 


53.2-43.1 




18 


7.9 










h;^co 


219.909 


3i,2- 2i,i 


0.41 


2.9 


1.0 


8.8 


0.84 


1.1 


1.6 


356.176 


52,3- 42,2 


<0.36 


2.4 


0.55 


6.3 


<0.41 


0.96 


0.92 




348.032 


52,3- 42,2 


0.15 


0.85 


<0.32 


2.6 


0.26 


<0.32 


<0.26 



- means frequency not observed. 



Table B.2. Observed Une fluxes / TuBdV (K km s"') for CH3OH in the 230 GHz window. 



Frequency 
[GHz] 


Transition 


AFGL 2591 


G24.78 


G75.78 


Sources 
NGC 6334 IRS 1 


NGC 7538 IRS 1 


W3(H20) 


W33A 


CH3OH 


218.440 


43,1,0- 32,1,0 




12.5 


3.0 










219.994 


23l4.9.0- 22i4.8.0 


<0.36 


<0.49 


<0.26 


2.0 


<0.32 


<0.47 


<0.41 


220.079 


84.4.0- 74,3.0 


0.55 


4.1 


1.1 


18 


1.3 


3.6 


0.54 


227.815 


169,8.0.1- 159.7.0.1 




2.1 


0.60 




0.87 


1.8 




229.939 


19l2.7.0.2- 2O12.8.O.2 


<0.31 


0.56 


<0.32 


5.1 


0.47 


<0.38 


<0.25 


230.027 


3 1.3.0- 42,3.0 


0.53 


3.4 


0.68 


17 


1.6 


4.3 


0.92 


230.293 


222200- 2I3190 


<0.31 


0.38 


<0.32 


5.1" 


<0.41 


0.68 


<0.25 


230.369 


22l3.9.o- 21 13,8.0 


<0.31 


0.26 


<0.32 


4.0 


0.67 


0.22 


<0.25 


231.281 


106,4,0,2- 96,3,0,2 


0.34 


7.6* 


0.93 


13 


0.82 


4.5 


1.7 


232.417 


106,4,0,1- 96,3,0,1 


0.64 


2.8 


0.72 


10 


<0.38 




0.94 


232.784 


I811.8.O.I- 17ll.7.0.1 








17 








232.927 


14681- 14771 








1.7 








232.946 


IO470- 1 157O 








16 








233.097 


13671- 13/61 








0.83 








234.683 


43.1.0.2- 53.2.0.2 


0.62 


4.3 


1.3 


18 


1.4 


4.1 


2.2 


234.698 


5 1.5.0- 62.5.0.5 


<0.24 


<0.36 


0.31 


10 


0.64 


1.6 


1.0 


240.242 


54.1.0- 64,2.0.5 


0.50 


3.7 


0.85 


12 


0.79 


5.6 


0.92 


240.741 


263,23,0- 262,24,0 


<0.25 


0.66 


<0.3 


4.8 


0.47 


0.33 


<0.20 


240.960 


53,2,1,1- 43,2,1,1*^ 


0.54 


3.0^ 


0.61 


17 


1.1 


3.4 


0.94 


241.043 


226,16.0- 235.ig 


<0.25 


<0.40 


<0.3 


2.6 


0.45 


<0.54 


<0.26 


241.057 


93.6.0- lOo.io.o 


<0.25 


<0.40 


<0.3 


3.3 


<0.35 


<0.54 


<0.26 


248.282 


17l0.7.0,2- 17l0,8,0,l 










<0.35 


3.0 




248.885 


1610,7,0,2- 169,7,0,1 


0.60 


2.4 


0.87 




0.96 


2.7 


<0.20 


"CH3OH 


220.363 


173150- 180180 


<0.19 


0.15 


<0.36 


1.3 


<0.30 


<0.46 


<0.19 


223.332 


64,3,1,2- 75,3,1,2 


<0.25 






1.3 


<0.36 






225.571 


147,8,0- 136,8,0 


<0.32 


0.64 


<0.24 


0.84 


<0.32 


0.50 


<0.23 


239.169 


192170- 19oi90 




<0.49 


0.19 






<0.38 


<0.29 


241.036 


22o220- 2I318O 


<0.25 


<0.40 


<0.30 


0.83 


<0.35 


<0.36 


<0.42 



- means frequency not observed. 



" Blend withHCOOCH3. 
* Blend with CH3CN. 
^ Blend with 52,3,2- 42 3 2 • 
Blend with HCOOCH3. 
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Table B.3. Observed line fluxes J TMsdV (K km s"^) for CH3OH in the 345 GHz window and its isotopic species. 



Frequency 


Transition 








Sources 








[GHz] 




AFGL2591 


G24.78 


G75.78 


NGC 6334 IRS 1 


NGC 7538 IRS 1 


W3(H20) 


W33A 


CH3OH continued 


337.969 


/ /1 3 1 1 — d/i 111 


0.48 


3.2 


0.76 


14 


1.3 


3.6 


0.77 


338.125 


/ A A n — 0"^ 1 n 


2.9 


20 


4.7 


56 


5.9 


24 


5.8 


338.345 


' 3,4,U "3,4,U 


3.9 


35 


7.7 


71 


8.1 


30 


11 


338.409 


1 A A n i — 111 1 n 1 

'4,4,U,1 "3,J,U,1 


3.9 


39 


7.9 


59 




66 


11 


338.431 


/ 1 1 ci — fin ^ rt 

' 1 , /,U ^U,o,U 




4.0 


0.84 










338.442 






4.7 


1.4 










338.456 






5.6 


1.5 










338.475 


' 0,1, u ^0,1, U 




5.9 


1.8 










338.486 


/ 1 n 1 /T — 0^; 1 n 1 /9 

O.J,U,l/Z "0,l,U,i/Z 




6.9 


1.6 










338.504 


1 ~i A — Mien 
' z,o,U ^1,3, U 




6.8 


1.9 










338.513 


7^ 3 2— 6*? A t 
j,j,u,z 4,2,0,2 




14 


3.3 










338.530 


76,2,0- 65,1,0 




7.2 


1.8 










338.542 


75,2,0,1/2- 65,2,0,1 




17 


4.7 










338.560 


72,5,0- 62,5,0 




8.6 


2.1 










338.583 


75.2,0- 65,2,0 




10 


2.5 










338.615 


74,3,0- 64,3,0 




20 


8.0 










338.640 


75,3,0,1- 64,2,0,1 




13 


3.1 










338.722 


73/5,5/3,0- 62/4,4/2,0 




27 


6.5 










340.141 


22,0,0,1- 32,1,0,1 


0.89 


5.6 


1.6 


27 


3.2 


4.1 


1.7 


340.394 


1611,5,0,1/2- 17ll,6,0,l/2 


0.65 


2.2 


0.58 


19 


1.6 


1.1 


0.78 


355.964 


16io,7,i- 169,7,1 


<0.29 


<0.47 


<0.53 


5.6 


<0.41 


<0.43 


<0.68 


356.007 


158,7,0,2- 158,8,0,1 


2.4 


9.5 


3.9 


29 


6.8 


13 


3.4 


'^CHjOH 


345.084 


22,0,0,1- 32,1,0,1 


<0.23 


0.26 


<0.30 


4.1 


<0.33 


<0.29 


<0.19 


348.101 


116,6,0- 106,5,0 


<0.18 


2.0 


0.27 


10.0 


0.61 


<0.44 


0.19 



- means frequency not observed. 
" Blend with 76,2,o,i/2- 65,1,0,1/2. 



Table B.4. Observed Une fluxes / Tw&dV (K km s-^) for CH2CO. 



Frequency 


Transition 








Sources 








[GHz] 




AFGL 2591 


G24.78 


G75.78 


NGC 6334 IRS 1 


NGC 7538 IRS 1 


W3(H20) 


W33A 


81.586 


4l,3-3,,2 


0,32 




<0.75 




0.45 


0.16 




140.127 


7l,7- 61,6 










0.72 


0.85 




244.712 


12i,ii- 111, 10 


<0.39 




<0.78 




1.8 


1.6 




262.619 


132.12- 122,11 


<0.33 




<0.47 




<1.7 


1.5 




220.178 


lli,ii-10uo 


<0.19 


1.2 


0.41 


5.0 


0.63 


0.65 


0.37 


222.199" 


llo,ii- 10o,io 


<0.31 


1.3 


0.38 


5.9 


0.51 


0.87 


0.37 


222.229 


112,10- 102,9 


<0.31 


0.38 


<0.22 


2.7 


<0.41 


0.20 


0.11 


222.315 


1 12.9- 102,8 


<0.31 


0.20 


<0.22 


1.5 


<0.41 


<0.35 


<0.25 


240.187 


12l,12- 111, 11 


<0.24 


1.5 


0.50 


4.5 


0.46 


1.4 


0.51 


346.600 


17l.l6- 16i,l5 


<0.26 


1.4 


0.52 


8.4 


0.81 


1.4 


0.60 



- means frequency not observed. 

" Blend of lines Hon - lOoio, II39 - IO38, and II38- IO37. 



Table B.5. Observed Une fluxes / Tus.dV (K km s'^) for CH3CHO. 



Frequency 
[GHz] 


Transition 


AFGL 2591 


G24.78 G75.78 


Sources 
NGC 6334 IRS 1 


NGC 7538 IRS 1 


W3(H20) W33A 


98.901 


5l,4- 4l,3 


0.29 


<0.14 






0.22 


0.38 




61,6- 5l,5 










0.40 


0.50 


149.507 


81.8- 7l,7 


<0.15 


<0.36 






0.69 


0.54 


168.093 


9l.9- 81,8 


<0.21 


<0.30 






0.30 


0.48 



- means frequency not observed. 
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Table B.6. Observed Une fluxes / TuBdV (K km s'^) for C2H5OH. 



Frequency 
[GHz] 


Transition 


AFGL 2591 


G24.78 


G75.78 


Sources 
NGC 6334 IRS 1 


NGC 7538 IRSl 


W3(H20) 


W33A 


220.155 


243_22,2- 242,23,2 


<0.24 


<0.38 


<0.36 


1.1 


0.78 


<0.39 


<0.24 


222.138 


197.13/12,0— 187,11/12,1 


<0.39 


<0.43 


<0.32 


1.6 


<0.32 


<0.27 


<0.32 


225.279 


135,9,1— 125,8,1 


<0.41 


<0.43 


<0.24 


0.57 


<0.26 


<0.21 


<0.29 


225.283 


135,8,1- 125,7,1 


<0.41 


<0.43 


<0.24 


0.9 


<0.26 


<0.21 


<0.29 


225.457 


133,11,1- 123,10,1 


<0.41 


0.43 


<0.24 


0.66 


<0.26 


<0.21 


<0.29 


225.546 


134,10,1- 124,9,1 


<0.41 


<0.43 


<0.24 


0.56 


<0.26 


<0.21 


<0.29 


225.560 


133,11,0- 123,10,0 


<0.41 


<0.43 


<0.24 


1.9 


<0.26 


<0.21 


<0.29 


225.660 


134,10,0- 124,9,0 


<0.41 


<0.43 


<0.24 


0.85 


<0.26 


<0.21 


<0.29 


226.581 


262,24,2- 26i,25,2 


<0.30 


<0.48 


<0.21 


2.6 


<0.26 


<0.35 


<0.30 


230.954 


165,11,2- 164,12,2 


<0.32 


0.59 


<0.21 


2.9 


<0.29 


<0.32 


<0.42 


230.991 


14o,14,2- 13l,13,2 


<0.32 


0.19 


<0.21 


2.8 


0.34 


0.21 


<0.42 


232.035 


185,14,2- 184,15,2 


<0.30 


<0.43 


<0.26 


0.56 


<0.30 


_ 


<0.41 


232.319 


235,19,2- 234,20,2 


<0.30 


<0.43 


<0.26 


0.45 


<0.30 


_ 


<0.41 


232.405 


175,13,2- 174,14,2 


<0.30 


<0.43 


<0.26 


0.95 


<0.30 


_ 


<0.41 


232.491 


14o,14,0- 13o,13,0 


<0.30 


0.24 


<0.26 


0.99 


<0.30 


_ 


<0.41 


232.809 


165,12,2- 164,13,2 


_ 


_ 


<0.30 


2.0 


_ 


_ 


_ 


232.929 


145,9,2- 144,10,2 


_ 


_ 


<0.30 


0.54 


_ 


_ 


_ 


233.571 


135,8,2- 134,9.2 


<0.41 


<0.43 


<0.24 


0.87 


<0.25 


<0.21 


<0.29 


234.524 


255,21,2- 254,22,2 


<0.30 


<0.48 


<0.21 


1.6 


0.23 


<0.35 


<0.30 


234.666 


105,5,2- 104,6,2 


<0.30 


<0.48 


<0.21 


<0.62 


0.24 


<0.35 


<0.30 


234.715 


105.6.2- 104,7,2 


<0.30 


<0.48 


<0.21 


2.0 


0.14 


<0.35 


<0.30 


234.853 


95,4,2- 94,5,2 


<0.30 


<0.48 


<0.21 


1.5 


<0.25 


<0.35 


<0.30 


234.984 


85,3,2- 84,4,2 


<0.30 


<0.48 


<0.21 


0.84 


<0.25 


<0.35 


<0.30 


234.992 


85,4,2- 84,5,2 


<0.30 


0.38 


<0.21 


1.9 


<0.25 


<0.35 


<0.30 


240.110 


272,25,2- 27l,26,2 


<0.30 


<0.36 


<0.26 


<0.72 


<0.30 


0.16 


<0.41 


240.782 


152,13,0- 14l,13,l 


<0.32 


<0.53 


<0.30 


1.1 


<0.28 


<0.45 


<0.26 


240.839 


14l,13,0- 13l,13,l 


<0.32 


<0.53 


<0.30 


2.0 


<0.28 


<0.45 


<0.26 


338.099 


187,11,2- 186,12,2 


<0.30 


1.7 


<0.20 


2.8 


<0.25 


<0.30 


0.14 


338.110 


187,12,2- 186,13,2 


<0.30 


<0.29 


<0.20 


2.2 


<0.25 


<0.30 


0.15 


338.163 


102,8,|- 9l,8,0 


<0.30 


0.37 


<0.20 


1.1 


<0.25 


<0.30 


<0.29 


339.511 


97,3/2,2- 96,4/3,2 


<0.30 


0.64 


<0.21 


1.3 


<0.28 


<0.39 


<0.21 


339.544 


87,2/1,2- 86,3/2,2 


<0.30 


<0.26 


<0.21 


2.8 


0.32 


<0.39 


<0.21 


339.566 


77,0/1,2- 76,1/2,2 


<0.30 


<0.26 


<0.21 


1.6 


<0.28 


<0.39 


<0.21 


339.979 


94.6.2- 83,5.2 


<0.30 


1.2 


<0.21 


6.1 


0.73 


<0.39 


<0.21 


340.189 


65.1/2.2- 54.2/1.2 


<0.36 


3.4 


0.68 


15 


1.6 


1.7 


1.1 


340.420 


94.5.2- 83,6.2 


<0.36 


<0.55 


<0.26 


4.0 


0.38 


<0.36 


0.78 


345.173 


77.0/1.0- 66.0/1.1 


<0.29 


<0.43 


<0.30 


2.4 


<0.26 


<0.21 


<0.24 


345.295 


21 121.1- 2O1.20.1 


<0.29 


0.55 


<0.30 


2.5 


<0.26 


<0.21 


<0.24 


345.408 


21 0.2 1.1- 2O0.20.1 


<0.29 


0.30 


<0.30 


5.1 


<0.26 


0.51 


<0.24 


346.620 


208,12/13,0- 198,11/12,0 


<0.26 


0.10 


<0.36 


1.3 


0.34 


<0.32 


<0.23 


346.817 


207,13/14,0- 197,12/13,0 


<0.26 


0.82 


<0.36 


3.9 


0.48 


0.22 


0.28 


346.929 


206,15,1- 196,14,1 


<0.26 


0.71 


<0.36 


1.4 


<0.19 


<0.32 


<0.23 


346.963 


21 0,2 1,2- 20 1,20,2 


<0.26 


2.0 


<0.36 


4.4" 


0.80 


0.89 


<0.23 


347.887 


205,15,0- 195,14,0 


<0.23 


0.38 


<0.32 


2.4 


0.23 


<0.32 


<0.26 


347.916 


204,17,0- 194,16,0 


<0.23 


2.1 


<0.32 


8.1 


1.4 


0.74 


<0.26 


347.975 


222,20,2- 213,19,2 


<0.23 


0.57 


<0.32 


2.7 


0.36 


<0.32 


<0.26 


348.720 


44,0/1,1- 33,0/1,0 


<0.36 


0.65 


<0.27 


<0.36 


<0.26 


<0.36 


<0.26 


348.847 


106,5/4,0- 95,5/4,1 


<0.36 


0.97 


<0.27 


<0.36 


<0.26 


<0.36 


<0.26 


352.858 


21 1,20,2- 202,19,2 


<0.32 


1.2 


<0.51 


4.9 


0.49 


<0.29 


<0.42 


353.035 


123,9,2- 112,10,2 


<0.32 


0.97 


<0.51 


5.2 


0.64 


0.58 


<0.42 


363.968 


217,15/14,1- 207,14/13,1 




0.50 


<0.38 










364.233 


217,15/14,0- 207,14/13,0 




0.57 


<0.38 











- means frequency not observed. 
"Blend with C2H5CN. 
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Table B.7. Observed Une fluxes / TuBdV (K km s'^) for HCOOH. 



Frequency 


Transition 








Sources 








[GHz] 




AFGL2591 


G24.78 


G75.78 


NGC 6334 IRS 1 


NGC 7538 IRS 1 


W3(H20) 


W33A 


223.916 


102.9-91.8 


<0.21 


— 


<0.17 


_ 


0.76 


0.19 


_ 


247.514 


1 15.6/7- 105.5/6 




_ 


_ 


— 


0.66 


2.4 


_ 


257.975 


12l,12- lll.ll 


<0.38 


_ 


<1.9 


_ 


<0.89 


<0.95 


_ 


262.103 


12o,i2- llo.ii 


- 


- 


- 


- 


<0.36 


<1.3 


- 


220.038 


lOo.io- 9o,9 


<0.19 


<0.29 


<0.29 


1.3 


0.37 


<0.46 


<0.19 


225.513 


103.7-93.6 


<0.32 


0.73 


<0.19 


2.6 


<0.25 


0.88 


0.40 


338.144 


154,12- 144,11 


0.23" 


2.1" 


<0.16 


3.7 


0.49" 


1.2« 


0.28 


338.202 


153,13- 143,12 


<0.24 


0.53 


0.26 


3.4 


0.74 


0.97 


0.24 


338.249 


154,1]- 144,10 


<0.24 


1.1 


<0.16 


3.1 


0.35 


0.72 


0.20 


345.031 


16o,16- 15o,15 


<0.23 


<0.32 


<0.20 


2.3 


0.51 


0.33 


<0.19 


346.719 


152,13- 142,12 


0.21 


0.49 


0.19 


2.9 


0.81 


0.36 


0.17 


356.137 


162,15- 152,14 


0.65 


0.28 


0.34 


1.8 


0.83 


0.55 


<0.54 



- means frequency not observed. 
"Blend with C2H5CN. 



Table B.8. Observed Une fluxes / TuBdV (K km s'^) for HNCO. 



Frequency 


Transition 








Sources 








[GHz] 




AFGL 2591 


G24.78 


G75.78 


NGC 6334 IRS 1 


NGC 7538 IRS 1 


W3(H20) 


W33A 


HNCO 


219.656 


103,7/8- 93,6/7 


<0.29 


<0.41 


<0.31 


1.2 


<0.25 


<0.47 


<0.41 


219.734 


102,9-92,8 


<0.29 


0.59 


0.77 


3.2 


0.72 


0.90 


0.78 


219.737 


102,8-92,7 


<0.29 


1.1 


a 


a 


a 


a 


a 


219.798 


10o,io-9o,9 


0.88 


4.4 


1.7 


8.4 


0.65 


1.4 


3.0 


240.876 


lli,ii-10i,io 


<0.25 


2.7 


1.3 


4.3 


0.55 


0.99 


1.4 


352.898 


161,15- 15l,14 


0.89 


3.3 


2.0 


7.7 


1.5 


2.0 


2.1 










HN'^ 


CO 








219.803 


IO2.9- 92.8 


<0.36 


<0.41 


<0.26 


<0.98 


<0.32 


<0.39 


<0.41 


240.881 


lli,ii-10i,io 


<0.25 


<0.46 


<0.36 


<0.41 


<0.28 


<0.54 


<0.26 



- means frequency not observed. 



" Blend with HNCO line at 219734 MHz. 



Table B.9. Observed Une fluxes / TuhdV (K km s'^) for NH2CHO. 



Frequency 
[GHz] 


Transition 


AFGL 2591 


G24.78 


G75.78 


Sources 
NGC 6334 IRS 1 


NGC 7538 IRS 1 


W3(H20) 


W33A 


227.606 


1 lo.u- lOo.io 




0.90 


<0.20 




a 


<0.39 




232.274 


1 12,10- 102,9 


<0.5 


<0.27 


0.15 


2.2 


0.32 




<0.60 


233.529 


1 16,5/6- 1064/5 


<0.15 


<0.32 


<0.19 


0.73 


0.56 


<0.21 


<0.42 


233.595 


1 15,6/7- 105,5/6 


<0.15 


0.64 


<0.19 


i> 


0.93 


1.4 


0.56 


233.736 


114,8- 104,7 


<0.15 


0.56 


<0.19 


1.3 


0.23 


<0.21 


0.30 


339.716 


168,8/9- 158,7/8 


<0.3 


1.3 


<0.17 


3.8 


1.2 


0.22 


<0.34 


339.780 


167,9/10- 157,8/9 


<0.3 


0.93 


<0.17 


2.2 


0.95 


<0.39 


<0.34 


339.904 


166,10/11- 156,9/10 


<0.3 


1.2 


<0.17 


4.1 


0.81 


<0.39 


<0.34 


340.491 


163,14- 153,13 


<0.36 


1.3 


0.18 


3.0 


0.70 


0.45 


0.75 


340.536 


164,13- 154,12 


<0.36 


0.55 


0.31 


2.3 


0.49 


<0.36 


<0.36 


345.183 


17o,17- 16o,16 


<0.29 


1.2 


<0.24 


2.1 


0.66 


0.44 


0.34 


345.327 


16l,15- 15l,l4 


<0.29 


<0.32 


<0.24 


3.8 


2.2 


<0.21 


<0.38 



- means frequency not observed. 
"Blend with C2H5 OH. 

* Blend with C2H5OH at 225404 MHz in other sideband. 
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Table B.IO. Observed line fluxes / TMsdV (K km s-^) for CH3CN. 



Frequency 


Transition 








Sources 








[GHz] 




AFGL 2591 


G24.78 


G75.78 


NGC 6334 IRS 1 


NGC 7538 IRS 1 


W3(H20) 


W33A 


CH3CN 


220.476 


128.0,K- 1 Is.CL 


<0.24 


0.45 


<0.43 


1.0 


<0.54 


<0.53 


<0.24 


238.913 


137,0,K- 127,0,L 




1.4 


<0.25 






<0.50 


0.38 


238.972 


136,0,K- 126,0,K 




2.8 


0.47 






1.8 


0.73 


239.064 


134,0,K- 124,0,L 




3.8 


0.83 






1.9 


0.97 


239.097 


133,0,K- 123,0,L 




5.9 


1.2 






3.8 


1.5 


239.120 


132.0,K" 122,0,L 




5.0 


1.1 






3.4 


1.5 


239.133 


13i,o.K- 12o/l,O.L 




5.7 


1.8 






8.1 


1.4 


239.138 


13o.().K^ 12()_(),L 




5.2 


0.98 






<0.50 


0.69 


240.090 


13l,2,K- 12l,2,L 


<0.3 


0.94 


<0.31 


3.0 


a 


0.29 


<0.32 


348.645 


19ll,0,K- 18ll,0,L 


<0.36 


<0.36 


<0.43 


0.64 


<0.59 


<0.58 


<0.24 


CH^'CN 


220.485 


126,0,K- ll6,0,L 


<0.24 


0.30 


<0.43 


0.35 


<0.54 


<0.53 


<0.24 


238.905 


135,O.K- 125,0,L 




0.17 


<0.25 






<0.50 


<0.34 


238.946 


134,0,K- 124,0,L 




0.57 


<0.25 






<0.50 


<0.34 


238.978 


133,0,K- 123,0,L 




1.06 


<0.25 






<0.50 


0.38 


239.001 


132,0.K- 122,0,L 




0.73 


<0.25 






<0.50 


<0.34 


348.853 


198,().K- ISs.O.L 


<0.36 


0.36 


<0.43 


<0.36 


<0.59 


<0.58 


<0.24 



- means frequency not observed. 

K refers to main transition of the upper level i.e. 12, 13, or 19 +1, +0, or -1 

L refers to the main transition of the lower level i.e. 11, 12 , or 18 +1, +0, or -1. 

"Blend with CH3OH. 



Table B.ll. Observed line fluxes f TuiidV (K km s"') for C2H5CN. 



Frequency 
[GHz] 


Transition 


AFGL 2591 


G24.78 


G75.78 


Sources 
NGC 6334 IRSl 


NGC 75 3 8 IRS 1 


W3(H20) 


W33A 


218.390 


24321- 233,20 




0.71 


<0.41 










219.506 


242.22- 232,21 


<0.36 


0.85 


<0.26 


2.2 


<0.32 


<0.39 


<0.41 


222.918 


25 1,24- 24i,23 


<0.32 






3.4 


<0.36 






223.385 


26l,26- 25l,25 


<0.32 






1.9 


<0.36 






227.781 


253,22- 243,21 




1.8 


<0.26 




<0.32 


<0.39 




228.483 


25223- 242,22 


<0.24 


1.4 


<0.36 


2.4 


<0.30 


<0.33 


<0.24 


231.310 


26i 25— 25i 24 


<0.32 


0.79 


<0.21 


2.5 


<0.36 


0.61 


<0.42 


233.443 


265,22- 255,21 


<0.15 


1.2 


<0.24 


1.0 


<0.41 


0.94 


<0.29 


233.498 


265,21- 255,20 


<0.15 


1.5 


<0.24 


1.9" 


<0.41 


0.63 


<0.29 


240.319 


28l,28- 27i,27 


<0.30 


1.4 


<0.26 


3.2 


<0.38 


0.66 


<0.41 


240.429 


28o,28- 27o,27 


<0.30 


1.1 


<0.26 


<0.58 


<0.38 


<0.47 


<0.41 


248.781 


29l,29- 28l,28 


<0.32 


1.1 


<0.30 




<0.35 


0.60 


<0.26 


248.869 


29o,29- 28o,28 


<0.32 


1.2 


<0.30 




<0.35 


1.1 


<0.26 


339.895 


392,38- 382,37 


<0.30 


1.2 


<0.21 


2.8 


<0.36 


<0.39 


<0.21 


339.968 


382,36- 372,35 


<0.30 


1.4 


<0.21 


2.9 


<0.36 


<0.39 


<0.21 


340.149 


39l,38 - 38l,37 


<0.36 


1.0 


<0.26 


<0.62 


<0.41 


<0.36 


<0.23 


340.433 


38i2,26/27- 37i2,25/26 


<0.36 


0.63 


<0.26 


1.1 


<0.41 


<0.36 


<0.23 


340.483 


38l0,28/29- 37(0,27/28 


<0.36 


0.83 


<0.26 


1.4 


<0.41 


<0.36 


<0.23 


340.551 


38l5,23/24- 37i5,22/23 


<0.36 


<0.35 


<0.26 


0.63 


<0.41 


<0.36 


<0.23 


340.576 


389,29/30- 379,28/29 


<0.36 


0.85 


<0.26 


<0.62 


<0.41 


<0.36 


<0.23 


346.712 


258,17/18- 257,18/19 


<0.26 


<0.47 


<0.36 


1.9 


<0.24 


<0.32 


<0.23 


346.864 


228,14/15- 227,15/16 


<0.26 


<0.47 


<0.36 


0.96 


<0.24 


<0.32 


<0.23 


346.874 


384,34- 374,33 


<0.26 


0.82 


<0.36 


2.4 


<0.24 


<0.32 


<0.23 


346.898 


218,13/14- 217,14/15 


<0.26 


0.49 


<0.36 


1.8 


<0.24 


<0.32 


<0.23 


346.925 


208,12/13- 207,13/14 


<0.26 


<0.47 


<0.36 


0.93 


<0.24 


<0.32 


<0.23 


348.261 


392,37- 382,36 


<0.23 


1.8 


<0.21 


3.6 


<0.27 


<0.32 


<0.26 


348.345 


40i,39- 39i.38 


<0.23 


1.2 


<0.21 


3.5 


<0.27 


<0.32 


<0.26 


353.089 


233,20- 222,21 


<0.32 


0.27 


<0.51 


0.58 


<0.30 


<0.29 


<0.42 


353.235 


1 16,5/6- 105,6/5 


<0.32 


<0.23 


<0.51 


2.9 


<0.30 


<0.29 


<0.42 



- means frequency not observed. 
" Blend with NH2CHO. 
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Table B.12. Observed Une fluxes / TMsdV (K km s-^) for CH3CCH. 



Frequency 


Transition 








Sources 








[GHz] 




AFGL 2591 


G24.78 


G75.78 


NGC 6334 IRS 1 


NGC 7538 IRS 1 


W3(H20) 


W33A 


222.014 


136- 126 


<0.23 


<0.27 


<0.25 


0.51 


<0.41 


<0.27 


<0.25 


222.061 


135- 125 


<0.23 


0.18 


<0.25 


<0.59 


<0.41 


<0.27 


<0.25 


222.099 


134- 124 


<0.23 


0.44 


<0.25 


1.2 


0.19 


0.49 


<0.25 


222.129 


133- 123 


0.63 


2.5 


0.54 


7.6 


0.92 


2.0 


0.58 


222.150 


132- 122 


0.87 


2.2 


0.81 


7.2 


1.2 


2.1 


0.89 


222.163 


13,-12i 


1.2 


4.3 


1.0 


8.1 


1.7 


1.7 


1.9 


222.167 


13o- 12o 


1.5 


3.4 


1.5 


8.1 


1.2 


3.2 


2.0 


239.088 


146- 136 




0.46 


<0.17 






<0.32 


<0.34 


239.138° 


145- 135 




0.69 


0.98 






<0.32 


<0.34 


239.179 


144- 134 




0.62 


<0.17 






0.70 


<0.34 


239.211 


143- 133 




2.6 


1.4 






1.9 


0.86 


239.234 


142- 132 




3.2 


1.4 






8.0* 


1.3 


239.248 


14i- 13i 




4.0 


1.6 






1.5 


1.9 


239.252 


14o- 13o 




4.4 


2.1 






3.4 


2.0 



- means frequency not observed. 
" Blend with CH3CN. 

* Blend with OCS at 231062 MHz in the main sideband. 



S. E. Bisschop et al.: Testing grain-surface chemistry in massive hot-core regions 
Table B.13. Observed line fluxes / T^BdV (K km s-^) for HCOOCH3 for the 230 GHz band. 



Frequency 
[GHz] 


Transition 


AFGL 2591 


G24.78 


G75.78 


Sources 
NGC 6334 IRS 1 


NGC 7538 IRS 1 


W3(H20) 


W33A 


211.537 


176,11- 166,10 E 


<0.36 


_ 


— 


8.6 


_ 


_ 


<0.14 


211.575 


176,11- 166,10 A 


<0.36 


_ 


_ 


6.3 


_ 


_ 


0.19 


211.771 


175,13- 165,12 E 


<0.36 


_ 


_ 


6.2 


_ 


_ 


<0.14 


211.785 


175,13- 165,12 A 


<0.36 


_ 


_ 


6.2 


_ 


_ 


<0.14 


218.281 


17i,i6- 16i,i5 E 


_ 


0.48 


<0.41 


_ 


_ 


_ 


_ 


218.297 


173,14- 163,13 A 


_ 


0.46 


<0.41 


_ 


_ 


_ 


_ 


220.167 


172,15- 164,12 E 


<0.24 


0.82 


<0.36 


7.1 


0.42 


0.98 


<0.24 


220.190 


174,13- 164,12 A 


<0.24 


0.67 


<0.36 


6.4 


0.66 


<0.33 


<0.24 


221.979 


189,9/10- 179,8/9 A 


<0.39 


0.71 


<0.32 


5.7 


0.84 


1.1 


<0.32 


223.125 


187,11- 177,10 E 


<0.32 


_ 


_ 


7.2 


0.30 


_ 


_ 


223.135 


187,12- 177,1 1 


<0.32 


_ 


_ 


8.9 


<0.43 


_ 


_ 


223.163 


187,11- 177,10 A 


<0.32 


_ 


_ 


<1.3 


0.14 


_ 


_ 


225.609 


193,16- 183,15 E 


<0.41 


<0.27 


<0.24 


4.9 


<0.32 


1.2 


0.39 


225.619 


193,17- 183,16 A 


<0.41 


1.0 


<0.24 


5.7 


<0.32 


1.9 


0.18 


226.713 


202,18- 192,18 E 


<0.30 


0.65 


<0.21 


9.1 


0.56 


0.60 


<0.30 


226.719 


202,19- 192,18 A 


<0.30 


1.3 


<0.21 


14 


0.66 


0.73 


0.52 


226.773 


203,18- 193,17 E 


<0.30 


<0.30 


<0.21 


<1.4 


0.44 


<0.35 


<0.30 


226.779 


20u9- 19,,i8 A 


<0.30 


<0.30 


<0.21 


<1.4 


0.53 


<0.35 


<0.30 


227.561 


21 1,21- 20i,2o E" 


_ 


5.3 


0.73 


_ 


1.4 


6.1 


0.34 


233.247 


19l6,3/4- 18l6,2/3 A 


<0.41 


<0.27 


<0.24 


1.5 


<0.32 


<0.21 


<0.29 


233.256 


19l6,3- I816.2 E 


<0.41 


<0.27 


<0.24 


3.0 


<0.32 


<0.21 


<0.29 


233.269 


19l6,4- 18l6,3 E 


<0.41 


<0.27 


<0.24 


1.3 


<0.32 


<0.21 


<0.29 


233.310 


19l5,4/5- 18l5,3/4 A 


<0.41 


<0.27 


<0.24 


3.1 


<0.32 


<0.21 


<0.29 


233.316 


19l5,4- 18l5,3 E 


<0.41 


<0.27 


<0.24 


2.0 


<0.32 


<0.21 


<0.29 


233.331 


19l5,5- 18l5,4 E 


<0.41 


<0.27 


<0.24 


2.3 


<0.32 


<0.21 


<0.29 


233.395 


19l4,5/6- 18l4,4/5 A* 


<0.41 


<0.27 


<0.24 


4.0 


<0.32 


<0.21 


<0.29 


233.414 


19l4,6- 18l4,5 E 


<0.41 


<0.27 


<0.24 


2.1 


<0.32 


<0.21 


<0.29 


233.505 


19,3,6- 18l3,5 E^ 


<0.41 


1.4 


<0.24 


6.1 


<0.32 


1.6 


<0.29 


233.525 


19,3,7- 18,3,6 E 


<0.41 


0.72 


<0.24 


2.8 


<0.32 


<0.21 


0.40 


233.649 


19,2,7- 18|2,6 E 


<0.41 


0.22 


<0.24 


<0.74 


<0.32 


<0.21 


<0.29 


233.655 


19,2,7/8- 18l2,6/7 A 


<0.41 


1.5'' 


<0.24 


4.4 


<0.32 


<0.21 


<0.29 


233.671 


19l2,8- 18l2,7 E 


<0.41 


0.18 


<0.24 


2.3 


<0.32 


<0.21 


0.17 


233.754 


182,16- 172,15 E 


<0.41 


1 .1 


<0.24 


7.5 


<0.32 


<0.21 


<0.29 


234.502 


199,10/11- 189,9/10 A 


<0.30 


0.94 


<0.21 


7.3 


0.41 


<0.35 


<0.30 


240.021 


19i,,8-18,,,7E 


<0.30 


1.1 


<0.26 


6.8 


0.72 


<0.42 


<0.41 


240.035 


193.,6- 183.,5 A 


<0.30 


1.1 


<0.26 


7.0 


0.68 


1.1 


0.29 


248.617 


207,,4- 197,,3 A 


<0.32 


0.69 


<0.30 




<0.35 


0.91 


<0.26 


248.634 


207,,4- 197,,3 E 


<0.32 


0.94 


<0.30 




<0.35 


0.96 


<0.26 


248.744 


207,,3- 197,,2 E 


<0.32 


<0.33 


<0.30 




0.69 


<0.45 


<0.26 


248.750 


232,22- 22i,22 E*" 


<0.32 


3.8 


<0.30 




1.3 


2.0 


<0.26 


248.787 


207,13- 197,12 A 


<0.32 


<0.33 


<0.30 




0.24 


<0.45 


<0.26 


249.031 


205,16- 195,15 E 


<0.32 


0.51 


<0.30 




<0.35 


0.33 


<0.26 


249.047 


205,16- 195,,5 A 


<0.32 


0.77 


<0.30 




<0.35 


0.59 


<0.26 



- means frequency not observed. 

" Blend with 2l2,2o- 202,,9 E,21 , 20, .20 A, and 21o,2i- 20o,2o A. 

* Blend with line at 233397 19i4,5- 18i4,4 E. 

' Blend with 19,3 6- 18,3,5 A, and 19i3,7- 18i3,6 A. 

Blend with C2H5CN. 
" Blend with 23i,23- 22i22 E, 232,22- 222,21 E, 23i,23- 222,21 E, 



40 S. E. Bisschop et al.: Testing grain-surface chemistry in massive hot-core regions 

Table B.14. Observed Une fluxes / TMsdV (K km s-^) for HCOOCH3 in the 345 GHz band. 



Frequency 
[GHz] 


Transition 


AFGL 2591 


G24.78 


G75.78 


Sources 
NGC 6334 IRS 1 


NGC 7538 IRS 1 


W3(H20) 


W33A 


338.356 


278,19- 268,18 E 


<0.30 


<0.30 


<0.20 


8.8 


0.23 


1.6 


<0.29 


338.397 


277,21- 267,20 E 


<0.30 


<0.30 


0.13 


11 


0.85 


<0.30 


<0.29 


340.221 


31 13,19- 31 12,19 E" 


<0.36 


<0.35 


<0.26 


2.7 


<0.41 


<0.36 


<0.23 


345.069 


28i4,i4— 27 14,13 E* 


<0.29 


1.7 


<0.30 


14 


<0.33 


1.5 


<0.24 


345.090 


28i4,i5- 27i4,i4 E 


<0.29 


0.69 


<0.30 


8.2 


<0.33 


0.78 


<0.24 


345.281 


18l3,5/6- 18l2,6/7 A 


<0.29 


<0.27 


<0.30 


3.5 


<0.33 


<0.21 


<0.24 


345.462 


28i3,i5- 27i3,i4 E 


<0.29 


<0.27 


<0.30 


4.8 


<0.33 


<0.21 


<0.24 


345.467 


28i3,i5/i6- 27)3, 14/15 A 


<0.29 


1.1 


<0.30 


8.9 


<0.33 


0.32 


<0.24 


345.486 


28i3,i6- 27i3,i5 E 


<0.29 


0.16 


<0.30 


5.7 


<0.33 


<0.21 


<0.24 


345.975 


28i2,i6- 27i2,i5 E 


_ 


_ 


_ 


_ 


0.68 


<0.29 


_ 


345.985 


28l2,16/17- 27i2,15/16 A 


_ 


_ 


_ 


_ 


0.78 


2.0 


_ 


346.001 


2812,17- 27i2,i6 E 


— 


_ 


_ 


— 


0.39 


0.64 


— 


346.660 


2811,17- 27ii,i6 E 


<0.26 


<0.47 


<0.36 


7.5 


0.67 


<0.32 


<0.23 


346.675 


28ii, 17/18- 27n, 16/17 A 


<0.26 


1.1 


<0.36 


11 


0.67 


<0.32 


<0.23 


346.687 


28ii,i8- 27ii,i7 E 


<0.26 


0.27 


<0.36 


6.7 


0.31 


<0.32 


<0.23 


348.050 


284,24- 274,23 E 


<0.23 


1.8 


0.11 


12 


0.73 


1.5 


<0.26 


348.066 


286,23 - 276,22 A 


<0.23 


1.6 


<0.32 


11 


0.60 


1.3 


<0.26 


348.909 


289,20— 279 19 E 


<0.23 


1.2 


<0.32 


<0.43 


<0.33 


<0.36 


<0.30 


348.915 


289,20- 279,19 A 


<0.23 


1.5 


<0.32 


7.7 


<0.33 


<0.36 


<0.30 


356.215 


2921,8/9- 2821.7/8 A 


<0.36 


<0.47 


<0.53 


2.6 


<0.41 


<0.36 


<0.68 


356.297 


2920,9/10- 2820,8/9 A 


<0.36 


0.49 


<0.53 


6.3 


0.30 


<0.36 


<0.68 


356.404 


29i9,io/ii- 2819,9/10 A 


<0.36 


<0.47 


<0.53 


4.2 


<0.41 


<0.36 


<0.68 



- means frequency not observed. 

" Blend with 31 13,15- 31 12,19 A, and 3 113,19- 31 12,20 A. 

* Blend with 281444- 271443 A, and 281445- 271444 A. 



Table B.15. Observed line fluxes / T„bdV (K km s'^) for CH3OCH3. 



Frequency 
[GHz] 


Transition 


AFGL 2591 


G24.78 


G75.78 


Sources 
NGC 6334 IRSl 


NGC 7538 IRSl 


W3(H20) 


W33A 


222.029 


212,20,2/3- 21l,21,2/3 


<0.39 


0.26 


<0.32 


3.1 


<0.41 


0.42 


<0.25 


222.035 


212,20,0/1- 21l,21,0/l 


<0.39 


0.72 


<0.32 


2.5 


<0.41 


0.37 


<0.25 


222.239 


43,2,2- 32,1,2 


<0.39 


0.16 


<0.32 


1.6 


<0.41 


<0.27 


<0.25 


222.248 


43,2,1/3- 32,1,1/3 


<0.39 


0.86 


<0.32 


5.3 


<0.41 


0.37 


0.37 


222.255 


43.2,0- 32,1,0 


<0.39 


0.53 


<0.32 


4.7 


<0.41 


0.50 


0.18 


222.259 


43,1,1/2- 32,1,1/2 


<0.39 


0.60 


<0.32 


2.6 


<0.41 


0.48 


0.18 


222.326 


25 3,23,a:- 244,20.0 


<0.39 


0.15 


<0.32 


2.3 


<0.41 


<0.27 


<0.25 


225.599 


12i,i2,a:- llo.ii.s: 


<0.41 


3.1 


0.73 


11.8 


0.83 


4.0 


1.0 


226.347 


14i,i3,a:- 132,i2,a: 




1.5 


<0.41 










227.648 


265,21,*:- 264,22,*: 




0.40 


<0.26 




<0.32 


<0.39 




227.658 


243,22,0/1- 242,23,0/1 




0.81 


<0.26 




<0.32 


<0.39 




230.140 


172,15,*:- 163,14,/r 


<0.39 


1.7 


<0.32 


3.3 


<0.41 


1.4 


<0.25 


233.633 


255,20,*:- 254,21,A- 


<0.15 


1.1 


,0.24 


2.6 


<0.32 


0.46 


<0.23 


240.978 


53,3,2- 42,2,2 


<0.32 


a 


<0.30 


2.6 


<0.35 


<0.45 


<0.26 


240.985 


53,3,1- 42,2,1 


<0.32 


3.4° 


<0.30 


8.3 


0.41 


2.8 


0.56 


240.990 


53,3,0- 42,2,0 


<0.32 


a 


<0.30 


3.3 


<0.35 


<0.45 


0.19 



- means frequency not observed. 
K is either 0,1,2, or 3. 

" 53,3,1- 42,2,1 transition is blended with 53,3,2- 42,2,2 and 53,3,0- 42,2,o- 



